Astronomy & Astrophysics manuscript no. HFI' Spectral' Bands 


©ESO 2013 


March 22, 2013 





Planck 2013 results. IX. HFI spectral response 

Planck Collaboration: R A. R. Ade^^, N. Aghanim^^, C. Armitage-Caplan^^ M. Arnaud^^, M. Ashdown^^'^, R Atrio-Barandela^^ J. Aumont^"^, 

C. Baccigalupi^^ A. J. Banday^^'^^, R. B. Barreiro^^, E. Battaner^^ K. Benabed^^'^^ A. Benoit^^ A. Benoit-Levy^^'^^'^^ J.-R Bernard^^, 
M. Bersanelli32'45, R Bielewicz86'lO'^^ J. Bobin6^ J. J. Bock^^'^^ J. R. Bond^ J. Borrill^^'^^ E R. Bouchet^^'^^ E Boulanger5\ M. Bridges63,6,58^ 
M. Bucher\ C. Burigana"^^'^^, J.-E Cardoso^^'^'^^ A. Catalano^^'^^ A. Challinor^^'^^'^^ A. Chamballu^^'^^'^'^, R.-R. Chary^\ X. Chen^\ 
L.-Y Chiang^^ H. C. Chiang^^'^, R R. Christensen^^'^^ S. Church^^^ D. L. Clements^^ S. Colombi^^'^^ L. R L. Colombo^^ ^i , c. Combet^^ 
B. Comis6^ E Couchot6\ A. Coulais^^ B. R Crill^i'^^ A. Curto^'^^ E Cuttaia4\ L. Danese^^ R. D. Davies^^ R de Bernardis^i, A. de Rosa4\ 
G. de Zotti4l'^^ J. Delabrouille\ J.-M. Delouis^^'^^ E-X. Desert4^ C. Dickinson^^^ J. M. Diego^^ H. Do\q^^^^\ S. Donzelli^^ O. Dore^^'^^ 
M. Douspis^\ X. Dupac^^ G. Efstathiou^^ T. A. EnBlin^\ H. K. Eriksen^^ E. Ealgarone^^ E Einelli4^'4^ O. Eorni^^'^^ M. Erailis^^ 
E. Eranceschi44, S. Galeotta^^ K. Ganga\ M. Giard^^'^^ Y. Giraud-Heraud\ J. Gonzalez-Nuevo^^js^ ^ M. Gorski^^'^^^ S. Gratton63,58^ 
A. Gregorio^^'^^, A. Gruppuso"^^, E K. Hansen^^, D. Hanson^^'^^'^, D. Harrison^^'^^, S. Henrot-Versille^"^, C. Hernandez-Monteagudo^^'^\ 
(T^ D. Herranz^^, S. R. Hildebrandt^\ E. Hivon^^'^^ M. Hobson^ W. A. Holmes^\ A. Hornstrup^^ W. Hovest^\ K. M. Huffenberger^^ 

^ G. Hurier54'6^ T. R. Jaffe^^'^^ A. H. Jaffe^^ W. C Jones2^ M. Juvela^^, E. Keihanen^^, R. Keskitalo^^ i^ T. S. Kisner^^ R. KneissP^ 8, 

O J- Knoche^\ L. Knox^^, M. Kunz^^'^^^ H. Kurki-Suonio2'^'3^ G. Lagache^"^, J.-M. Lamarre^^ A. Lasenby^'^^ R. J. Laureijs^^ C. R. Lawrence^\ 
CN J. R Leahy^^ R. Leonardi^^ C. Leroy^'^'^^'^^, J. Lesgourgues^^'^^ M. Liguori2^ R B. Lilje^^ M. Linden-V0rnle^^ M. Lopez-Caniego^^, 

^ R M. Lubin27, J. R Mac^as-Perez6^ B. Maffei^^^ N. Mandolesi^^'^'^^ M. Maris^^ D. J. Marshall6^ R G. Martin^ E. Martinez-Gonzalez^^ 
^ S. Masi^\ S. Matarrese^'^, E. Matthai^\ R. Mazzotta^^, R. McGehee^\ A. Melchiorri^^ "^^, L. Mendes^^, A. Mennella^^'^^ M. Migliaccio^^'^^ 
^ S. Mitra49'6i, M.-A. Miville-Deschenes^^^ A. Moneti^^ L. Montier^^'^^ G. Morgante4\ D. Mordock^^ D. Munshi^^^ J. A. Murphy^^ 

^ R Naselsky^4,35^ p ^^^^1^ p Natoli^o^^^ c. B. Netterfieldl^ H. U. N0rgaard-Nielsenl^ C. North^^ E Noviello^^ D. Novikov^^ I. Novikov^\ 
O S. Osborne^^ C. A. Oxborrow^^ E. Raci^^ L. Ragano^^'^^ E. Rajot^"^, D. Raoletti'^'^''^^ E. Rasian'^^ G. Ratanchon\ O. Rerdereau^"^, L. Rerotto^^ 
^ E Perrotta^^ E Piacentini3\ M. Piat\ E. Pierpaoli^^, D. Pietrobon^^ S. Plaszczynski^^, E. Pointecouteau^^'^^ G. Polenta^'^^ N. Ponthieu^4'4^ 

, I L. Popa5^ T. Poutanen39'24,2 q ^^ Prsiii^^, G. Prezeau^^'^^ S. Prunet^^'^^ J.-L. Fuget^\ J. R Rachen^O'^^, M. Reinecke^\ M. Remazeilles^^ ^ 

^ C. Renault^^ S. Ricciardi"^^, T. Riller^\ I. Ristorcelli^^'^^, G. Rocha^^'^\ C. Rosset\ G. Roudier^'^^'^\ B. Rusholme^\ D. Santos^^ G. Savini^^ 
E. R S. Shellard^^ L. D. Spencer*^^^ J.-L. Starck6^ V. Stolyarov^'^^'^^ R. Stompor\ R. Sudiwala^^ E Sureau6^ D. Sutton^^'^^ 
A.-S. Suur-Uski24'39^ j._F. Sygnet^^ J. A. Tauber^^, D. Tavagnacco43,33^ p. Terenzi^^, M. Tomasi^^ M. Tristram6\ M. Tuccii^'6\ G. Umana^^ 
4^ L. Valenziano4\ J. Valiviita3'^'24'^^ B. Van Tent6^ R Vielva^^ E Villa4\ N. Vittorio3\ L. A. Wade^^ B. D. Wandelt^^'^^'28, D. Yvon^^ 

JpH A. Zacchei^^ and A. Zonca^^ 

O (Affiliations can be found after the references) 

C/^ Received XX, 2013; accepted XX, 2013 



ABSTRACT 



> 

o 

o 
in 



The Planck High Erequency Instrument (HEI) spectral response was determined through a series of ground based tests conducted with the HEI 
focal plane in a cryogenic environment prior to launch. The main goal of the spectral transmission tests is to measure the relative spectral response 
(including the level of out-of-band signal rejection) of all HEI detectors to a known source of electromagnetic radiation individually. This was 
determined by measuring the interferometric output of all detection channels for radiation propagated through a continuously scanned polarizing 
Eourier transform spectrometer. As there is no on-board spectrometer within HEI, the ground-based spectral response experiments provide the 
CO definitive data set for the relative spectral calibration of the HEI. Knowledge of the relative variations in the spectral response between HEI 
detectors allows for a more thorough analysis of the HEI data. The spectral response of the HEI is used in Planck component separation and 
CO data analysis, this includes extraction of CO emission observed within Planck bands, dust emission, Sunyaev Zeldovich sources, and intensity to 
^— H polarization leakage. The HEI spectral response data have also been used to provide unit conversion and colour correction analysis tools. 

While previous papers have already described the pre-flight experiments conducted on the Planck HEI, this paper focuses on the analysis of the 
, ^ pre-flight spectral response measurements and the derivation of data products, e.g. band-average spectra, unit conversion coefficients, and colour 
k>j correction coefficients, all with related uncertainties. Verifications of the HEI spectral response data are provided through comparisons with HEI 
ffight data. This validation includes use of HEI Zodiacal emission observations to demonstrate out-of-band spectral signal rejection better than 
^ 10^. The accuracy of the HEI relative spectral response data is verified through comparison of complementary bandpass based and flight-data 
based unit conversion coefficients and colour correction coefficients. These coefficients include those based upon HEI observations of CO, dust, 
and Sunyaev-Zeldovich emission. General agreement is observed between the ground-based spectral characterization of HEI and corresponding 
in-flight observations, within the quoted uncertainty of each; explanations are provided for any discrepancies. 

Key words. Astronomical instrumentation, methods and techniques - Instrumentation: detectors - Instrumentation: photometers - Space vehicles: 
instruments - Cosmology: observations - cosmic background radiation 



1. Introduction 

This paper, one of a set associated with the 2013 release of 

data from the Planck^ mission (Planck Collaboration I 2013), 

* Corresp. author: L.D.Spencer, Locke . Spencer@astro .cf.ac.uk 

^ Planck (http://www.esa.int/Planck) is a project of the 
European Space Agency (ESA) with instruments provided by two sci- 
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describes the determination and verification of the Planck High 
Frequency Instrument (HFI) spectral response. As the HFI em- 
ploys a series of broad-band photometric receivers, an accurate 
understanding of the relative spectral response of each detector 
within a frequency channel, and that of each frequency chan- 
nel within the instrument, is important in data processing and 
analysis (Planck Collaboration VI 2013); this is particularly im- 
portant for Planck component separation (Planck Collaboration 
XII 2013) where the magnitude of foreground components is 
much greater than the uncertainties associated with the cosmic 
microwave background (CMB) signals. 

The Planck High Frequency Instrument (HFI) spectral re- 
sponse was determined through a series of ground based tests 
conducted with the HFI focal plane in a cryogenic environment 
prior to launch. One of the main goals in pre-flight calibration 
testing was to measure the relative spectral response (includ- 
ing the level of out-of-band signal rejection) of all HFI detec- 
tors to a known source of electromagnetic (EM) radiation indi- 
vidually. This was determined by measuring the interferometric 
output of all detection channels for radiation propagated through 
a continuously scanned polarizing Fourier transform spectrome- 
ter (FTS). As there is no on-board spectrometer within HFI, the 
ground-based spectral response experiments provide the defini- 
tive spectral calibration of HFI. Pre-flight component level spec- 
tral characterization testing is described in detail in Ade et al. 
(2010). The pre-flight system level spectral characterization test- 
ing is described in detail in Pajot et al. (2010). This paper will 
discuss the testing itself as needed to provide context, but will 
primarily concentrate on the spectral characterization measure- 
ments, their analysis, and their utility within the HFI consortium 
and the Planck legacy data archive. Additional details of the HFI 
pre-flight calibration are found in Rosset et al. 2010, Lamarre 
et al. 2010, and Mafl'ei et al. (2010). The average spectrum for 
each of the HFI bands is illustrated in Figure 1 (see Sect. 3.1 for 
details). 

An accurate understanding of the spectral response of HFI 
is critical in HFI data processing and analysis. Information de- 
rived from the spectral calibration is important in many aspects 
of component separation (Planck Collaboration XII 2013). Unit 
conversion factors and colour corrections, which are important 
when dealing with signals of varying spectral profiles within a 
photometric channel, are derived for HFI using the transmis- 
sion spectra. Spectral mismatch between detectors within a given 
band must be understood in order to accurately interpret multi- 
detector averages and diff'erences that otherwise yield systematic 
errors and increased uncertainty in data products. This is espe- 
cially important in the evaluation of weak components including 
polarized signal. 

This paper presents the propagation of the raw spectral 
characterization data through its processing and analysis to 
yield detector level spectral response data in Sect. 2. Sect. 3 
then describes the analysis of the detector spectral transmission 
data to provide advanced spectral data products such as band- 
average spectra and unit conversion / colour correction coeffi- 
cients. Sect. 4 presents the evaluation of the HFI spectral re- 
sponse and data products through comparisons with flight data. 
This includes comparisons with HFI observations of Zodiacal 
light (Planck Collaboration XIV 2013), CO emission (Planck 
Collaboration XIII 2013), Sunyaev-Zeldovich sources (Planck 



entific consortia funded by ESA member states (in particular the lead 
countries France and Italy), with contributions from NASA(USA) and 
telescope reflectors provided by a collaboration between ESA and a sci- 
entific consortium led and funded by Denmark. 
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Fig. 1. Band-average spectral transmission for each of the HFI 
frequency channels. 



Collaboration XXVIII 2013, Planck Collaboration XXIX 2013, 
and Planck Collaboration XXI 2013), and dust emission (Planck 
CoUaboration 2013, Planck CoUaboration XVII 2011, and 
Planck CoUaboration XIX 2011). 

Details of the spectral response of the Planck LFI instrument 
may be found in Planck Collaboration II (2013) and in Zonca 
et al. (2009). Further details on technical aspects of the HFI spec- 
tral response data are also available in the Planck Explanatory 
Supplement (Planck Collaboration ES 2013). 

2. Measurements 

This section outlines the pre-flight data collected allowing the 
HFI spectral response determination. The processing of the 
recorded data is outlined, and the resultant spectra are presented. 

2.1. Spectral Response Data Collection 

The Planck HFI detector spectral response data were measured 
using a broadband mercury arc lamp radiation source, a polar- 
izing FTS, an integrating sphere, and a rotating filter wheel, all 
coupled with the HFI focal plane in its evacuated cryostat. A 
reference bolometer (at approx. 4 K) was mounted within the in- 
tegrating sphere as an external measure of the radiation incident 
on HFI (see Fig. 3, Pajot et al. 2010). With this configuration, 
the entire HFI focal plane and the reference bolometer were ex- 
posed to the FTS modulated spectral signal synchronously. The 
ratio of the HFI detector spectra with the reference bolometer 
spectra provides a relative spectral response with the system- 
atics due to the test configuration removed. The HFI spectral 
response measurements were collected during June and July of 
2006 at the Institut d'Astrophysique Spatiale (IAS) laboratories 
in Orsay, France. For each HFI detector, roughly 100 interfero- 
grams were recorded with consistent scan parameters to allow 
a spectral resolution of approximately 0.5 GHz and a Nyquist 
frequency (Nyquist (1928) and Shannon (1948)) of approxi- 
mately 12000 GHz in the resultant spectra. The rotating filter 
wheel was placed before the entrance to the integrating sphere 
with two settings used in these observations (Pajot et al. 2010). 
A 10 cm" ^(approx. 300 GHz) ^ low-pass edge (LPE) filter was 
used for observations of the 100-217 GHz detectors while a 
36 cm" ^(approx. 11 000 GHz) LPE filter was used for the 353 - 

^ The term wavenumber will be used to refer to units of cm~^ 
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857 GHz observations. The LPE filtering external to HFI and the 
integrating sphere allowed better performance of the reference 
bolometer over the low frequency channel spectral range, as will 
be shown below. Further details of the experimental setup, in- 
cluding diagrams of the FTS, integrating sphere, and HFI focal 
plane locations, are found in Pajot et al. (2010). Additional rele- 
vant information is also found in Ade et al. (2010). 

There were two significant additional tests in the deriva- 
tion of the HFI detector spectral response beyond the scope 
of the IAS HFI FTS measurements. Optical efficiency exper- 
iments, also conducted at IAS, provide optical efficiency esti- 
mates for each HFI detector. When coupled with the FTS spec- 
tra this allows an estimate of the absolute spectral transmis- 
sion. More details on the optical efficiency tests are in Catalano 
(2008); Catalano et al. (2006). Additional filter measurements 
were recorded in the AIG test facility at Cardiff' during filter 
stack production (Ade et al. 2010). These measurements extend 
the IAS FTS spectral measurements far beyond the HFI spectral 
passband up to approximately 20 THz. The independent measure 
of the filter stack transmission is used for a portion of the HFI 
detector spectral transmission for regions of the spectrum where 
it is deemed to be of better quality then the IAS FTS measure- 
ments (i.e., for frequencies outside of the HFI band edge filter 
cut-off). The fflter stacks for each of the frequency bands are 
comprised of 5 ff Iters. There is an additional low frequency cut- 
on fflter for the 545 and 857 GHz bands as the waveguide cut-on 
is too low for the desired multi-mode performance (Ade et al. 
2010 and Murphy et al. 2010). 

2.2. Spectral Response Data Processing 

The raw detector signals were combined with a bolometer model 
(Holmes et al. 2008, Lamarre et al. 2010, and Planck HFI Core 
Team 2011) to both convert the signal into physical units and 
perform a detector nonlinearity correction (Nay lor et al. 2009). 
The recorded interferograms were processed and Fourier trans- 
formed individually (Bell 1972 and Davis et al. 2001), includ- 
ing phase correction (Forman et al. 1966 and Brault 1987) and 
apodization (Naylor & Tabic 2007). The resultant spectra were 
then averaged together to provide a mean and standard deviation 
for every independent spectral data point. Similar analysis was 
conducted for the reference bolometer measurements. The ra- 
tio of each detector average spectrum against the corresponding 
reference bolometer average spectrum was taken to obtain the 
relative spectral transmission. Fig. 2 demonstrates an example 
set of raw and average spectra for the HFI 100 GHz la detector. 
The signal-to-noise (S/N) ratios for the average spectra for each 
HFI detector and the relevant spectral region of the reference 
bolometer spectra are shown in Fig. 3. In terms of the noise of the 
HFI detector relative transmission spectra, the reference bolome- 
ter is the limiting case, especially for the 100 GHz detectors. It 
is important to note, however, that each HFI detector spectrum 
within a given frequency band is divided by the same reference 
bolometer spectrum. Therefore, the relative uncertainty between 
HFI detectors is indicated by the HFI-only marks of Fig. 3, even 
though the absolute uncertainty is dictated by the limiting S/N of 
the reference bolometer. The detector relative transmission spec- 
tra are normalized to have a maximum value of unity, with the 
optical efficiency test results (Catalano 2008) providing a mul- 
tiplicative term to obtain the absolute spectral response, i.e., the 
product of the normalized spectral response and the optical ef- 
ffciency factor provides an estimate of the absolute spectral re- 
sponse of a given detector. As the reference bolometer accepts 
In sr. of incident radiation within the integrating sphere, the rela- 




Frequency [GHz] 

Fig. 2. Sample bolometer spectra from HFI 100 GHz detector 
la. The black data represent individual spectra while the colored 
data represent the average of the approx. 100 individual spectra. 
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Fig. 3. HFI detector spectral S/N for the scan-averaged HFI de- 
tector spectra prior to taking the ratio with the scan-averaged 
reference spectra. Also shown is the S/N for the scan-averaged 
reference bolometer spectra over the same spectral region. 

live transmission spectrum for each HFI detector is also through- 
put ^ normalized by virtue of the reference bolometer ratio. 

Once the detector and reference spectra ratios have been ob- 
tained, the out-of-band spectral regions are modiffed to improve 
the overall data quality. A waveguide model is used for the low 
frequency pre-cut-on region and the component level fflter stack 
spectra are used for the post-cut-off' spectral regions (see Ade 
et al. (2010) for more details). For the 545 and 857 GHz chan- 
nels, the component level fflter spectra are also used in conjunc- 
tion with the waveguide model for the pre-cut-on spectral re- 
gion. Fig. 4 illustrates an example composite spectrum showing 
the relative FTS spectrum, the waveguide model, the component 
level fflter spectra, and the corresponding spectral transition re- 
gions where the external spectra are spliced onto the in-band 
spectrum. 

At this stage, the frequency sampling for a given detector 
may be slightly diff'erent than those within the same frequency 
channel as each detector signal is processed independently. To 
allow easier intra-channel comparisons, all spectra are then in- 
terpolated onto a common frequency grid within each of the HFI 

^ Throughput, i.e., A Q, is defined as the area - solid angle product of 
a diffraction limited system (Born & Wolf 1999). 
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Table 1. HFI regions of CO rotational transmission 



Band 


CO transition 




over- sampled 


[GHz] 


(-^upper -^lower) 


[GHz] 


region [GHz] 



100 
Frequency | 



250 



100 

217 

353 

545 
545 

857 
857 
857 
857 



115.2712018 

230.5380000 

345.7959899 

461.0407682 
576.2679305 

691.4730763 
806.6518060 
921.7997000 
1036.9123930 



109.67 

219.34 

329.00 

438.64 
548.28 

657.89 
767.48 
877.04 
986.57 



115.39 

230.77 

346.15 

461.51 
576.85 

692.17 
807.46 
922.73 
1037.95 



Fig. 4. Example ratioed spectrum (black), filter spectrum (red), 
and waveguide model (blue) for bcOO - 100 GHz la. The tran- 
sition frequencies are shown by the vertical dashed lines. The 
out-of-band transmission is a conservative over-estimates as the 
ratioed spectrum is scaled by the optical efficiency prior to out- 
of-band grafting of the waveguide model and filter data. 

frequency channels. As all measurements within a frequency 
channel are conducted synchronously, and all detectors are refer- 
enced against the same reference bolometer, uncertainties intro- 
duced as a result of the common frequency interpolation are ex- 
pected to be negligible. The largest differences occur as a result 
of differences in band-edge location, and thus the transition fre- 
quency when the component level filter spectra (at a lower spec- 
tral resolution) are used in place of the IAS spectra. The common 
frequency sampling adopted attempts to mitigate the inclusion 
of additional spectral uncertainty by selecting the limiting case 
within a frequency band, e.g., the lowest filter- spectrum transi- 
tion point is selected to avoid the interpolation increasing the 
localized spectral resolution of any detector spectra. This may 
result in a minor degradation of the spectral resolution for some 
detectors, but does not result in presenting data at a higher spec- 
tral resolution than it was originally measured with."^ 

Regions of the HFI spectral range with sensitivity to CO 
emission are identified (Jet Propulsion Laboratory 2004). To as- 
sist with the use of the HFI spectral response data with HFI CO 
studies (see Planck Collaboration XIII (2013)), the signal in the 
HFI CO spectral regions (see Table 1) is interpolated by a factor 
of 10 using the FTS instrument line shape (ILS) (Connes 1961). 
Although the data are presented at higher resolution, the native 
spectral resolution of the underlying data does not increase with 
this interpolation, i.e., the number of independent data points 
is not increased and the FTS instrument line-width remains the 
same. While this interpolation provides a more accurate estimate 
of the spectral transmission within these regions than, e.g., lin- 
ear interpolation, it is important to note that these data do not 
correspond to independent sampled frequency data points. Thus, 
a flag column has been added to the spectral transmission pro- 
file data files to indicate whether a given data point originates 
from the actual data, or is a result of the ILS-based interpolation. 
The interpolation regions have been extended to include the fre- 
quency range of other CO isotopes (i.e., CO, ^^CO, C^^O, and 
C^^O) over a range of radial velocities (i.e., ±300kms"^). The 
original data points within the over-sampled region have been 
preserved, so a flag filter on the data will restore the independent 
data points easily. 

The CO interpolation, discussed below, is an exception to this. 



2.3. Spectral Response Data Products 

The HFI detector spectral response data products are available 
within the database instrument model (see Planck Collaboration 
ES (2013)) and within the Planck legacy archive (PLA). The 
data is comprised of the spectral frequency in units of both GHz 
and cm"\ the normalized spectral response (with its associated 
uncertainty), a CO interpolation flag, and meta-data including 
waveguide and filter transition regions, optical efficiency, and 
date/version information. The spectral normalization is such that 
the maximum value of any given spectrum is unity. An esti- 
mate of the absolute spectral transmission is obtained through 
the product of the optical efficiency parameter and the normal- 
ized spectrum for a given detector or frequency band (see Planck 
Collaboration ES (2013) for further details). Figs. 5-10 present 
the spectral response for each of the HFI detectors. Several diag- 
nostic parameters are determined for each HFI detector spectrum 
(and the band-average and detector sub-set spectra, see Sect. 3.1 
and Table 4). These parameters include the cut-on and cut-off 
frequency, the effective bandwidth, the central frequency, vari- 
ous effective frequencies, and the optical efliiciency, all of which 
are defined below. 

The following definitions are used in the evaluation of the 
spectral transmission profiles for the HFI detectors: 

- Cut-on, Von- The Cut-on frequency defines where the spec- 
tral band or the high-pass filter frequency dependence goes 
from the minimum to the maximum value. The cut-on po- 
sition of the filter is usually referred to as the position in 
frequency (or wavenumber) space where the smooth varying 
function in question reaches half of the maximum power. 

In the most general case, this definition encounters a problem 
when dealing with spectra that have small peaks and can os- 
cillate above and below the half-power mark in getting to the 
maximum in band power. IMoreover the maximum in band 
power can be a peak value, an average value, or the band 
itself can have a "sloped" shape that biases the actual defi- 
nition of half-maximum level. The cut-on frequency in this 
work is defined as the lowest frequency occurrence of half- 
maximum amplitude. 

- Cut-off, Vojf: Similar arguments apply to the cut-off as for the 
cut-on. The cut-off frequency in this work is defined as the 
highest frequency occurrence of half-maximum amplitude. 

- Bandwidth, Ay.- The FWHIM of the optical band, given by 

Av = Voff - Von . (1) 
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Central Frequency, Vcen- The central frequency is defined as 
the average of Von and Voff . A more useful parameter, how- 
ever, is the effective frequency , Veff. 

Effective frequency, Veffi Alternatively the eff'ective fre- 
quency is usually defined by weighing the spectra by the 
frequency itself and is analogous to the most-probable fre- 
quency and is given by 



Veff ■ 



J VT\v)dv 
f T\y)dy 



(2) 



where T\y) is the spectral transmission including the optical 
efficiency term s. 

Integrated optical efficiency, si^i'. The integrated optical ef- 
ficiency is obtained by integrating the spectral transmission 
across the entire measured frequency range, and dividing by 
the detector bandwidth as follows 



^Int 



T{v)dv 



Ay 



(3) 



where e is the relative optical efficiency discussed in Sect. 2.2 
above, 6:int is the integrated optical eflftciency introduced 
here, t{v) is the normalized spectral transmission, and Ay 
is the bandwidth as described in Eq. 1 above. 
Spectral index effective frequencies, Va,tf['- As in the eff'ective 
frequency case, this is equivalent to weighing the transmis- 
sion spectrum by the frequency while including a diff'erent 
spectral index, a. The eff'ective frequency including spectral 
index is given by 



where Vq is the nominal band reference frequency. • 



There is an important distinction between e and 6:int in that the 
former is a scaling term, which accompanies the normalized 
transmission spectra (and is meaningless on its own), and the lat- 
ter is intended to represent an eff'ective optical efficiency over the 
specified bandwidth, i.e., an equivalent spectral rectangle/tophat 
function. 

In addition to the individual detector spectra, several other 
data products have been prepared for distribution. These include 
band-average spectra (see Sect. 3.1, and also Table 4), unit con- 
version / colour correction coefficients (see Sect. 3.2), and unit 
conversion/colour correction, UcCC, software scripts to accom- 
pany HFI data (see Planck Collaboration ES (2013)). These data 
products are discussed in the following sections. Tables 2 and 
3 report the parameters derived from the HFI detector spectral 
transmission profile data products as defined above, for the HFI 
band-average spectra; similar results for all HFI detectors are 
available in Planck Collaboration ES (2013). The spectral in- 
dices chosen for Table 3 correspond to the IRAS SED (a = -1), 
a planetary SED {a = 2), and a dust SED {a = 4). 



^ This frequency is somewhat arbitrarily defined as it does not have 
to be equal to y^ff above, but is a matter of definition. The choice of 
Vc for Planck detectors and frequency channels is discussed further in 
Sect. 3.2.1. 
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Fig. 5. The detector spectral transmission profiles for the HFI 
100 GHz detectors. 



(4) 3. Results 

This section presents data products derived from the HFI de- 
tector spectral responses described above. This includes band- 
average spectra, along with unit conversion and colour correc- 
tion algorithms and coefficients. 



3.1. Frequency Channel-Average Spectra 

Frequency channel- average transmission spectra, i.e., band- 
average spectra, are derived to complement the HFI frequency 
channel maps of various components (Planck Collaboration XII 
2013). To produce band- average spectra, individual detector 
spectra of a given frequency channel are weighted by a de- 
tector scaling factor to mimic the proportional weighting ap- 
plied in the HFI mapmaking algorithms (Planck Collaboration 
VI (2013) and Planck Collaboration VIII (2013)). This scaling 
factor, i.e., w/, is based on relative noise levels, spectral response, 
and instrument scan strategy, all of which are described below. 
Although eff'orts were made to duplicate the mapmaking rou- 
tines, the determination of the band-average transmission spectra 
is similar, but not identical, to its mapmaking counterpart. The 
method used to determine the band-average transmission spectra 
is described below. The individual detector weights described in 
Sect. 3.1.1 are identical to those used in the mapmaking scripts. 
The divergence lies in the hit-map normalization (Eq. 6), and 
the CMB normalization (Sect. 3.1.3). The hit map normalization 
could fully reproduce the approach of the mapmaking routines 
if an average were produced for each map-pixel (see Sect. 3.4), 
rather than producing a single scalar w/ coefficient for each de- 
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Table 2. HFI spectral response diagnostic parameters 



Spectrum Von [GHz] y^ff [GHz] Ay [GHz] Vctn [GHz] y^ff [GHz] ej^t 



100-avg 


84, 


.4 


+ 


0.3 


117, 


.36 


+ 


0.05 


32, 


.9 


+ 


0.3 


100. 


,89 


+ 


0.14 


101, 


.31 


+ 


0.05 


0, 


.304 


+ 


0.003 


100-detsetl 


84, 


.77 


+ 


0.09 


117, 


.81 


+ 


0.05 


33, 


.03 


+ 


0.11 


101, 


,29 


+ 


0.05 


101, 


.43 


+ 


0.07 


0, 


.265 


+ 


0.002 


100-detset2 


84, 


.29 


+ 


0.18 


117, 


.14 


+ 


0.05 


32, 


.85 


+ 


0.19 


100, 


,72 


+ 


0.09 


101, 


.25 


+ 


0.06 


0, 


.321 


+ 


0.003 


143-avg 


119, 


.994 


+ 


0.018 


165, 


.76 


+ 


0.04 


45 


.76 


+ 


0.05 


142. 


,875 


+ 


0.020 


142, 


.709 


+ 


0.015 


0, 


.3669 


+ 


0.0006 


143-detsetl 


120, 


.05 


+ 


0.03 


160, 


.18 


+ 


0.09 


40 


.13 


+ 


0.10 


140. 


,12 


+ 


0.05 


141, 


.45 


+ 


0.03 


0, 


.4614 


+ 


0.0017 


143-detset2 


118, 


.95 


+ 


0.08 


164, 


.9 


+ 


0.8 


45 


.9 


+ 


0.8 


141. 


,9 


+ 


0.4 


142, 


.27 


+ 


0.02 


0, 


.379 


+ 


0.007 


143-SWBs 


120, 


.17 


+ 


0.03 


166, 


.308 


+ 


0.018 


46 


.14 


+ 


0.04 


143. 


,238 


+ 


0.018 


143, 


.96 


+ 


0.03 


0, 


.3123 


+ 


0.0007 


217-av2 


188, 


.892 


+ 


0.011 


253, 


.419 


+ 


0.007 


64 


.527 


+ 


0.013 


221. 


,156 


+ 


0.006 


221, 


.914 


+ 


0.005 


0, 


.33850 


+ 


0.00012 


217-detsetl 


183, 


.3 


+ 


0.3 


253, 


.606 


+ 


0.020 


70 


.3 


+ 


0.3 


218. 


,46 


+ 


0.13 


220, 


.548 


+ 


0.009 


0, 


.3053 


+ 


0.0011 


217-detset2 


182, 


.159 


+ 


0.013 


253, 


.592 


+ 


0.007 


71 


.433 


+ 


0.016 


217. 


,875 


+ 


0.007 


220, 


.614 


+ 


0.009 


0, 


.34838 


+ 


0.00018 


217-SWBs 


189, 


.02 


+ 


0.03 


253, 


.247 


+ 


0.013 


64 


.22 


+ 


0.04 


221. 


,136 


+ 


0.017 


222, 


.957 


+ 


0.008 


0, 


.3226 


+ 


0.0002 


353-avg 


306, 


.8 


+ 


0.6 


408, 


.22 


+ 


0.02 


101 


.4 


+ 


0.6 


357. 


,5 


+ 


0.3 


361, 


.289 


+ 


0.008 


0, 


.335 


+ 


0.002 


353-detsetl 


303, 


.582 


+ 


0.015 


406, 


.333 


+ 


0.017 


102 


.75 


+ 


0.02 


354. 


,957 


+ 


0.011 


359, 


.156 


+ 


0.011 


0, 


.29902 


+ 


0.00014 


353-detset2 


318, 


.885 


+ 


0.014 


407, 


.86 


+ 


0.02 


88 


.97 


+ 


0.03 


363. 


,372 


+ 


0.013 


360, 


.870 


+ 


0.012 


0, 


.28730 


+ 


0.00015 


353-SWBs 


306, 


.3 


+ 


0.4 


408, 


.81 


+ 


0.03 


102 


.5 


+ 


0.4 


357. 


,56 


+ 


0.18 


361, 


.921 


+ 


0.011 


0, 


.3575 


+ 


0.0013 


545-avg 


469, 


.5 


+ 


0.5 


640, 


.81 


+ 


0.03 


171 


.3 


+ 


0.5 


555. 


,2 


+ 


0.3 


557, 


.54 


+ 


0.03 


0, 


.2612 


+ 


0.0008 


545-detsetl 


466, 


.44 


+ 


0.02 


642, 


.36 


+ 


0.04 


175 


.91 


+ 


0.04 


554. 


,40 


+ 


0.02 


557, 


.86 


+ 


0.03 


0, 


.28031 


+ 


0.00013 


545-detset2 


470, 


.9 


+ 


0.3 


638, 


.52 


+ 


0.11 


167 


.6 


+ 


0.4 


554. 


,73 


+ 


0.17 


556, 


.85 


+ 


0.05 


0, 


.2143 


+ 


0.0005 


857-avg 


743, 


.9 


+ 


0.4 


989, 


.78 


+ 


0.08 


245 


.9 


+ 


0.4 


866. 


,8 


+ 


0.2 


862, 


.68 


+ 


0.05 


0, 


.2165 


+ 


0.0004 


857-detsetl 


736, 


.9 


+ 


0.7 


990, 


.38 


+ 


0.06 


253, 


.4 


+ 


0.7 


863. 


,7 


+ 


0.4 


863, 


.42 


+ 


0.06 


0, 


.2121 


+ 


0.0006 


857-detset2 


741, 


.79 


+ 


0.13 


987, 


.01 


+ 


0.10 


245, 


.22 


+ 


0.17 


864. 


,40 


+ 


0.08 


861, 


.74 


+ 


0.07 


0, 


.21419 


+ 


0.00017 



Table 3. HFI spectral response effective frequencies 



Spectrum 




=-1 


[GHz] 




.=2 [GHz] 




.=4 [GHz] 


100-avg 


100. 


,36 


+ 


0.05 


103, 


.24 


+ 


0.05 


105 


.25 


+ 


0.04 


100-detsetl 


100. 


,49 


+ 


0.07 


103, 


.35 


+ 


0.06 


105 


.34 


+ 


0.06 


100-detset2 


100. 


,31 


+ 


0.07 


103, 


.19 


+ 


0.06 


105 


.21 


+ 


0.05 


143-avg 


141. 


,362 


+ 


0.015 


145, 


.457 


+ 


0.014 


148 


.234 


+ 


0.013 


143-detsetl 


140. 


,11 


+ 


0.03 


144, 


.22 


+ 


0.02 


147 


.05 


+ 


0.02 


143-detset2 


140. 


,91 


+ 


0.02 


145, 


.05 


+ 


0.02 


147 


.90 


+ 


0.02 


143-SWBs 


142. 


,64 


+ 


0.03 


146, 


.63 


+ 


0.02 


149 


.28 


+ 


0.02 


217-avg 


220. 


,111 


+ 


0.005 


225, 


.517 


+ 


0.006 


229 


.096 


+ 


0.007 


217-detsetl 


218. 


,666 


+ 


0.009 


224, 


.312 


+ 


0.009 


228, 


.038 


+ 


0.010 


217-detset2 


218. 


,697 


+ 


0.009 


224, 


.429 


+ 


0.009 


228, 


.200 


+ 


0.010 


217-SWBs 


221. 


,241 


+ 


0.008 


226, 


.395 


+ 


0.008 


229 


.834 


+ 


0.010 


353-avg 


358. 


,563 


+ 


0.008 


366, 


.763 


+ 


0.009 


372, 


.192 


+ 


0.010 


353-detsetl 


356. 


,386 


+ 


0.011 


364, 


.744 


+ 


0.012 


370 


.302 


+ 


0.013 


353-detset2 


358. 


,409 


+ 


0.012 


365, 


.850 


+ 


0.012 


370 


.837 


+ 


0.013 


353-SWBs 


359. 


,158 


+ 


0.011 


367, 


.455 


+ 


0.012 


372, 


.930 


+ 


0.014 


545-avg 


552. 


,22 


+ 


0.05 


567, 


.596 


+ 


0.017 


576 


.778 


+ 


0.014 


545-detsetl 


552. 


,43 


+ 


0.06 


568, 


.12 


+ 


0.02 


577, 


.458 


+ 


0.017 


545-detset2 


551. 


,76 


+ 


0.08 


566, 


.48 


+ 


0.03 


575, 


.32 


+ 


0.02 


857-avg 


854. 


,69 


+ 


0.11 


877, 


.724 


+ 


0.018 


891 


.462 


+ 


0.016 


857-detsetl 


855. 


,33 


+ 


0.16 


878, 


.67 


+ 


0.02 


892 


.59 


+ 


0.02 


857-detset2 


853. 


,89 


+ 


0.17 


876, 


.53 


+ 


0.03 


890 


.03 


+ 


0.02 



3.1 .1 . Noise Weighting 

The HFI detector noise equivalent temperature (NET) estimates 
(Planck Collaboration VI 2013 and Planck Collaboration VIII 
2013) are used to weight the individual detector signal in av- 
eraging within data processing. An attempt at duplicating this 
behaviour is made to obtain multi-detector average spectra. To 
determine the relative weights of individual detectors within an 
average, the inverse square of the detector NETs is normalized 
such that the sum total within the desired detector grouping is 
unity as follows 

w.- = ^^^^^^ where W = ^(1/(NET)2) . (5) 

The detector NETs and the Wi factors can be found in Planck 
Collaboration ES (2013). Two detectors have been omitted from 
contributing towards the band-average spectra due to random 
telegraphic signal (RTS), i.e., popcorn noise: 143GHz-8 and 
545 GHz-3. The w/ factors introduced here represent a general 
weighting concept that will be expanded upon further. As a re- 
sult of this, factors strictly following Eq. 5 will be referred to as 
wnet / in subsequent sections. 

3.1 .2. Detector Channel-Map Contribution Weighting 

Furthermore, the NET can be scan-normalized using the 
individual-detector pixel-hit maps available as standard data 
products (these will be made publicly available in the final re- 
lease of Planck data if not earlier, further details on the hit-maps 
can be found in Planck Collaboration ES (2013)), i.e.. 



(Z^,0^m,/(^,0))/(NET,-)2 
^m,i = — , (6) 

tector. The CIMB normalization described here is analogous to ^ 

the dipole calibration done in the standard HFI mapmaking (but, where 1^ is a normalization term as described above, Hi(6, (p) 
again, not identical). represents the hit-map counts for a given detector, sky position. 
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Fig. 6. The detector spectral transmission profiles for the HFI Fig. 7. The detector spectral transmission profiles for the HFI 
143 GHz detectors. 217 GHz detectors. 



and a given map, m (e.g., full-survey, nominal- survey, survey 1, 
etc.), the Yue,(i) represents summing over the entire map. A 
similar approach could be taken where the Y^q^^ summation is 
omitted; instead of a single w/ factor for a given map, m, this 
would result in a map of weighting factors of the same spatial 
resolution as the map, i.e., w/(^, 0). An example of this, using 
the nominal survey and survey 1 results for the 857 GHz band, 
is shown in Fig. 11 ^ , where it is clear that the 857-4 detector 
contributes less to the maps. 

A global factor, wh net can be obtained from each Wi{6, (p) 
map by choosing the statistical mean, or median, of the map, 
by taking the peak of a histogram of the map values, or through 
choosing some other diagnostic method. Fig. 12 illustrates his- 
tograms based on the detector Wm^iiO, cp) maps for the HFI full 
survey, nominal survey, and, in the 100 GHz case, survey 1. 
Similar histograms were also computed for individual sky sur- 
veys (HFI completed 5 fully sky surveys), masked surveys 
(masking varying percentages of the galactic plane and bright 
sources), and detector sub-set maps (see Table 4); these have 
been omitted from Fig. 12 for clarity. For each HFI detector in- 
cluded in the plot, the dashed-dotted vertical lines (marked as 
NET) indicate the wnet / factors resultant from Eq. 5, while the 
long-dashed vertical bars (marked as H NET) indicate those re- 
sultant from Eq. 6; these coefficients are generated from the full 
survey w^jiO, (p) map histograms. Similar results are available 
for various subsets, all converging towards the full survey values 

^ This figure, and subsequent map figures within this paper, were 
produced using a modified version of the HEALPix software routines 
(Gorski et al. 2005). 



(more detailed figures for each of the HFI bands, and various 
data sub-sets, are found in Planck Collaboration ES (2013)). The 
full survey values are displayed as they present stronger con- 
vergence than any given subset of the data. The incorporation 
of the detector hit count into the band-average scaling factors 
ensures that the resultant frequency channel spectra best repre- 
sents the eff'ective transmission spectrum for a given map. This 
is demonstrated by the diff'erence in position of the NET and 
H NET markings on the figure. In other words, the Planck scan- 
ning strategy is an important consideration in determining the 
band-average spectra due to relative hit-counts and integration 
time changing for diff'erent detectors and sky positions. 

The 857 GHz example shown is a special case, as the 857- 
4 bolometer exhibits higher noise properties than it's counter- 
parts. The 857-4 detector contribution relative weight is thus 
much less than the other 857 GHz detectors. The 857-4 his- 
togram (Fig. 12.b) is not symmetric, with an inflated tail towards 
zero weight. As a result, the 857 Wmj scalar value is intentionally 
ofl'-peak by a small factor. 

The band-average spectra made available for distribution 
are thus based on the Wm,i values from the full survey maps. 
While the nominal band-average spectra are produced using the 
detector- weight histogram peak values, the efl'ects of using ofl'- 
peak histogram weights on the resultant band-average spectra 
were investigated. Using compatible weight factors (2 = 1) 
from the histogram tails produce band-average spectra that vary 
from the nominal band-average spectra at the percent level. 
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Fig. 8. The detector spectral transmission profiles for the HFI Fig. 9. The detector spectral transmission profiles for the HFI 
353 GHz detectors. 545 GHz detectors. 



3.1 .3. Photometric Bandpass Averaging 

A map from detector /, m/, is given by 



Mi = Ki 



1 



dv{A^\Ti{y)dI, [Kcmb], 



(7) 



where Ki is the photometric dipole calibration factor (e.g., 
KcmbW"^), the (1 + 77/)/2 fraction is used to distinguish be- 
tween spider-web bolometers (SWB)s and polarization sensitive 
bolometers (PSB)s, (AQ)y represents the telescope throughput, 
Tiiy) represents the normalized spectral transmission, 6/ repre- 
sents the optical efficiency, and dly represents the diff'erential 
source intensity (see Eq. 9). As riy) is throughput normalized, by 
virtue of the ratio of the HFI detector spectra against a reference 
bolometer spectrum within an integrating sphere (see Pajot et al. 
(2010)), the (AQ)vT(y) term is further reduced to {AQ)y^T{y), i.e., 
the product of the throughput at the nominal reference frequency 
and the throughput normalized transmission spectra. For SWBs, 
the (1 + 7]i)/2 coefficient is 1, and for PSBs it is 1/2. 

The temperature derivative of the Planck function, evaluated 
at the CMB temperature, i.e., the CMB dipole, is the principle 
calibration source for HFI. This is used for the 100- 353 GHz 
channels, while FIRAS data, accompanied by a planet-based re- 
normalization, is the calibrator used for the 545 and 857 GHz 
channels (Planck Collaboration VI 2013). Based on the well 
known Planck function, By(T, y), the CMB dipole signal is as- 



sumed to have the following form: 
dBy(T, y) 



K = 



dT 



r=2.7255 K 

2hv^ 



c2(exp \hvl{KT)\ - 1) 



hv \ 



exp [hy/(KT)] ' 
exp [hy/(KT)] - 1 ^ 



(8) 





\ ^ 1 


r=2.7255 K 


m^ srHzK 



and the diff'erential source intensity is thus given by 

W 



dly = b[,dTcMB 



m^ sr Hz 



(9) 



All of the terms outside of the integral in Eq. 7 can be replaced 
by a single constant. A/, as follows 



Ai = Ki^-^(A^)y,€i 



KcMBm^sr 



W 



(10) 



If the source dly is a dipole CMB spectrum, then the map itself 
should be given by m/ = dTcMB ^ therefore 



At = |JjyT^(y)/7;j = Kt^^(AQ)y,ei 



KcMBm^sr 



W 



. (11) 



The corresponding channel map, M, is given as a weighted com- 
bination of the detectors comprising that channel. For individual 
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Fig. 10. The detector spectral transmission profiles for the HFI 
857 GHz detectors. 




Nominal Survey 



Survey 1 



Fig. 11. Detector relative map contribution weight factors, i.e., 
Wm,i as above, for the 857 GHz detectors. The left half represent 
the nominal survey (a., c, e., g.), while the right half represent 
the survey 1 subset of the data (b., d., f., h.). Data corresponding 
to detectors 1 through 4 are grouped in rows, i.e., 857-1 is illus- 
trated in (a.,b.), and 857-2,-3, and -4 are found in (c.,d.), (e.,f.), 
and (g.,h.), respectively. The grey shaded regions indicate where 
there was no sky coverage. 



weightings of w/, the map, M, neglecting polarization eff'ects for 
now, is given by 



M=^^-^ [Kcmb]. 



(12) 



Expanding this expression to include Eqs. 7, 10, and 11, we get 



(13) 



div [Kcmb] 



It is clear that the transmission and source components of the 
spectrum are separable, with the eff'ective spectral response 
given by 



AT(y) : 



lliiWiAiTiiv)) 



KcMBm^sr 



W 



where A is an arbitrary scaling factor such that 
M = A JdyT(y)dI, [Kcmb]- 



(14) 



(15) 



Similar to the derivation of Eq. 11, for a CMB dipole source, 
M = dTcMB and 



KcMBm sr 



W 



(16) 



Equation 14 can be rearranged to solve for r(v) as follows 



T(V) 



f \ 

w/T/(y) 



(17) 



The above expression, however, contains the desired T(y) on both 
sides of the equation. Since the right hand instance of T(y) is 
within a frequency integral, and thus will only result in a scal- 
ing factor to be applied to the average transmission spectrum, 
an alternate strategy will be used. The CMB -normalized channel 
transmission spectrum, '^^mb^^^' defined as 



^CMB^^Norm 





( Y 




W,T,(V) 





(18) 



where Norm[/(x)] is defined as /(x)/max[/(x)]. A CMB-weight 
factor, wcMB is introduced to define the contribution of each 
detector to the bandpass average. The above w/ factors (Eqs. 5 
and 6) are coupled with the derivative of the CMB spectral func- 
tion to determine the noise/CMB -normalized scaling factors as 
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Fig. 12. Histograms of the Wm,i band-average spectra scaling fac- 
tors. Values are shown for 100 GHz (a.) and 857 GHz (b.) detec- 
tors, including the full, nominal, and individual surveys (some 
are omitted for clarity). The vertical bars represent the resultant 
weight factor both with and without the Planck sky coverage 
and/or hit-maps taken into consideration. 



follows 



(S,w,) \dvri{v)K 











max< 


i 


w;r,(v) 




(S,w,)(Jjvt,(v)Z.;|7-cmb) 


> 



. (19) 



In the above expression, w/ is a scalar factor unique for each HFI 
detector. This not need be the case, however, as this case may be 
more further generalized by allowing the detector weight factor 
Wi to vary across the sky as in Eq. 6. This generalization results 
in a wcMB iCd, (p) photometric weighting factor, i.e., the relative 
weights varies across the sky, and with relative integration time, 
etc. Fig. 12 illustrates histograms of the detector weight factors 
across the sky for the HFI full mission data, and various survey 
sub-sets of this data. The difference between the static Wi fac- 
tor and the histogram peak value is demonstrated by the vertical 
bars in the figure. Thus, the frequency band average transmission 
spectra are comprised of the individual detector spectra propor- 
tionately scaled for both the relative response to the CMB spec- 



trum, and the relative noise level within a given channel. The 
resultant band-average transmission spectra are shown in Fig. 1 
above. 

As the 545 and 857 GHz channels are calibrated using 
FIRAS data, and subsequently renormalized using planet ob- 
servations (Planck Collaboration VI (2013)), rather than using 
the CMB dipole directly, it is important to investigate the use of 
Eq. 19 in deriving the band-average spectra for these channels. 
A comparison using both the w/ and wcmb i scaling factors for 
the 545 and 857 GHz channels was thus conducted. While these 
diff'erences for the 100 - 353 GHz channels are at the level of 
a few percent, they are at the 0.3 - 0.5 % level for the 545 and 
857 GHz channels. Furthermore, the diff'erences do not exceed 
the respective uncertainty of the corresponding detector spectra. 
As this normalization removes the dependence of the individ- 
ual T/(y) values on an absolute calibration, i.e., the optical effi- 
ciency, this CMB -normalization in the band-average spectra is 
maintained for all of the HFI bands. 



3.2. Unit Conversion and Colour Correction 

This section presents the formula used to obtain unit conversion 
and colour correction coefficients for use with the HFI data, and 
the method used to derive the uncertainties on these coefficients, 
and the coefficients. 



3.2.1 . Unit Conversion and Colour Correction Philosophy 

Broad-band detection instruments, including photometric instru- 
ments using band-defining filters such as HFI, measure power 
collected via an instrument collecting area for unresolved/point 
sources, and power collected within a given throughput for 
extended sources. Although such an instrument directly mea- 
sures power absorbed by the detectors, it is convenient to relate 
this power measurement to either flux density (unresolved/point 
sources) or specific intensity/brightness (extended sources) such 
that the combined spectral and throughput integrated signal is 
equal to the measured power. Thus, observation data expressed 
in units of brightness, specific intensity, or flux density, are in- 
trinsically associated with an assumed reference frequency and 
spectral energy distribution (SED) profile. 

Spectral calibration of broad-band photometric instruments 
is performed by observation of a source of known SED. Provided 
that observed sources have a similar SED (within the spectral 
band) to that of the calibration source, measurements are cali- 
brated by the ratio of the two observations. The general case is 
that observed sources have a diff'erent SED to that of the cal- 
ibration source(s). Any instrument observation is then related 
to a calibration observation by expression of both in terms of 
an equivalent intensity at a specified reference frequency. The 
equivalent intensity is defined by knowledge of the source SED 
for a given observation. A colour correction (Griffin et al. 2013) 
is used to relate measurements of one SED to those of another. 
There are two equivalent approaches to astronomical colour cor- 
rection. One approach is to determine the eff'ective frequency 
that corresponds to the assumed SED and measured intensity, 
and determine a diff'erent reference frequency, based upon a dif- 
ferent SED, for any other SEDs of interest. In converting be- 
tween SED types, the intensity remains the same, but at a dif- 
ferent reference frequency. The other approach is to determine 
the relative intensity for a given reference frequency, so the ref- 
erence frequency remains the same for various SEDs, but the 
intensity will vary. 
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The Planck HFI uses two calibration schemes (Planck 
Collaboration VIII 2013), one based on the differential CMB 
dipole spectrum, and another based upon the more local astro- 
physical sources. The 100-353 GHz channels are calibrated on 
the CMB dipole, which follows a dByldT]^^^^ SED profile, 

where data are provided in units of diff'erential CMB tempera- 
ture, i.e., KcMB- The 545 and 857 GHz channels are calibrated 
on a combination of galactic emission (typically dust following 
an approx. y"^ SED profile) and planetary emission (of an ap- 
prox. SED profile), where data are provided in units of bright- 
ness/intensity, i.e. MJysr"^ Furthermore, the 545 and 857 GHz 
data are scaled to equate with an SED following the IRAS con- 
vention (Beichman et al. 1988), which has a SED profile of the 
form y/y = constant. It is thus important to express observation 
data in both formats for various aspects of data analysis. This in- 
volves both unit conversion, where data are presented in a diff'er- 
ent unit, but remain consistent with a given SED (e.g., MJy sr"^ 
can be expressed as an equivalent brightness in K), and colour 
correction, where data are expressed with respect to a diff'erent 
assumed SED at the same reference frequency (e.g., changing 
from KcMB to MJy sr"^ with a SED of spectral index -1). 

In the millimetre - sub millimetre region of the electro- 
magnetic spectrum, the diff'use emission is made of several 
components whose summation represents the observed signal. 
Component separation algorithms (Planck Collaboration XII 
2013), when separating physical components with diff'erent SED 
profiles (e.g., emission from the CMB, thermal dust, spinning 
dust, free-free sources, synchrotron sources, CO lines, etc.), use 
models or templates for these components. The component sep- 
aration then must resolve an inverse problem going from sev- 
eral maps of broad-band measurements at diff'erent frequencies 
to component maps. The model is adjusted to minimize residu- 
als through fitting the sum of the components to the measured 
intensity. It is therefore impractical to have the various com- 
ponents all expressed at diff'erent eff'ective frequencies. Thus, 
Planck adopts a fixed reference frequency unit conversion/colour 
correction where the intensity is corrected for the assumed (or 
measured) SED of the source or component. 

Using KcMB calibration for the submillimeter channels, 
especially 857 GHz, should be avoided. This is because the 
submillimetre IRAS to CMB unit conversion depends heav- 
ily on the low-frequency region of the bandpass spectrum (see 
Fig. 14), which is known with less confidence than the main 
band. Therefore, the conversion of the 100-353 GHz data from 
a CMB to IRAS SED is less error-prone. 

Details on the derivation of the Planck unit conversion and 
colour correction coefficients are provided in the following sec- 
tion. 



3.2.2. Coefficient Formula Derivation 

The following conversion factors are derived for the individual 
HFI detectors and the frequency-channel average spectra: 

1. Convert [MJy/sr] (IRAS) ^ [KcmbJ. 

2. Convert [MJy/sr] (IRAS) ^ [Krj]. 

3. Convert [yszl ^ [KcmbI- 

4. Colour correction (power-law spectra and modified black- 
body spectra). 

5. CO correction. 

In general, the unit conversion terms are arrived at by equat- 
ing changes in intensity, expressed in various forms. Starting 



with the general expression J/y = J/y; e.g., 

J/y = (dIJdXiXdXi) = (dIJdXjXdXj) 



(20) 



Each side may also be multiplied by the spectral transmission, 
T(y), and integrated across the frequency band as follows 

/.vr(v)(^).Z, = /.vr(v)(^).X, . (21) 

This can be simplified to the form of a unit conversion coefficient 



dX, 
dX 



(22) 



In converting to/from Kcmb . the derivative of the Planck func- 
tion at the CMB monopole temperature (T^cmb = 2.7255 K, see 
Fixsen (2009)) will be used (see Eq. 8 and 9). 

Conversion to MJy/sr (IRAS) is accomplished using the 
IRAS convention (Beichman et al. 1988) ^ , v J/y = constant, 
such that 



W 



m^ sr Hz 



(23) 



where Vq is the band centre frequency and die is the effective 
intensity at the central frequency for a source of spectral index 
minus one. 

As Rayleigh-Jeans brightness temperature units are of- 
ten more convenient (i.e., human friendly) than that of 
Wm"^ sr"^ Hz"^ or even MJy sr"\ a ffux density to brightness 
temperature unit conversion is provided using the following re- 
lation (Rybicki & Lightman 1986): 



dTh 



Iv^k 



dly [Kb] 



(24) 



where K^, is the temperature expression of flux density, i.e., 
brightness temperature does not imply a Rayleigh-Jeans spectral 
profile. The brightness temperature unit of is selected over 
Krj to avoid any confusion between this definition, and that of a 
source exhibiting a Rayleigh-Jeans SED profile. 

For SZ conversion the following intensity expression, based 
on the Kompaneets non-relativistic SZ formula (Kompaneets 
(1957), Sunyaev & Zeldovich (1980), and Grainger (2001)), is 
used 



dly SZ = {K)(T) 



hv \ exp {hvl{kT)\ + 1 



kT j exp [hvlikT)] - 1 



-4 



Cvsz) 



^CMB 



W 



m^ sr Hz 



(25) 

where, again, b'^ is the temperature derivative of the Planck func- 
tion. 

In power-law colour correction, it is assumed that the source 
intensity follows a power law over the spectral region of interest, 
i.e., dly oc y^, providing the following expression 



dly a 



die a 



W 



m^ sr Hz 



(26) 



http : //lambda . gs£c . nasa . gov/product/iras/docs/exp . 
sup/ch6/C3.html 
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A colour correction to a modified-blackbody of the form 
dly oc v^By(y, T) (By defined above) is given by 



dL 



By(V, T) 



By(y,,T) 



dL 



w 



m sr Hz 



(27) 



For molecular rotational transitions, such as the CO Jl-0, 
J9-8 transitions, the desired specific intensity term is an eff'ec- 
tive brightness temperature, /STqq, in units of Kkms"^ For a 
Doppler line profile, v is equal to yco(l + which is closely 

approximated by vcoCl - vie) fox v « c. The intensity can be 
said to be distributed across a narrow velocity distribution, dv ^ 
, such that the integral over all velocities yields the temperature- 
velocity eff'ective brightness. To relate the eff'ective brightness 
across a narrow frequency range, dv, to frequency units rather 
than those of velocity, the relationship dvjdv = -vqoIc is used, 

Arco(^) 

Arcolv= [K] . (28) 

dv 

It is important to note that the eff'ective brightness temperature, 
i.e., AT^co^ is given in units of K km s"^ while this brightness 
distributed across a defined velocity interval {dv), or a defined 
frequency range {dv), i.e., Arcoly. is given in units of K. The 
CO transition intensity is thus given by the following relation 

dico = (ATcoDb'^j 

(29) 

= (ATco)r-^K \^^] ■ 

\cdv/ ^ Lm^srHz. 
The Rayleigh-Jeans(R-J) approximation temperature derivative 
used above is given as follows 

_ dBy^{T,v) 
^RJ ~ 



dT 



2v^k 



RJ 



w 



m sr Hz Krj 



(30) 



As the CO transitions occur at discrete frequencies, with line 
widths much narrower than the spectral resolution of the detector 
spectral transmission profiles, the CO intensity integral can be 
approximated by a delta function at the CO frequency, i.e., 

r ivT(v)dico = T(vco) (Arco) (— ) b'^l 

J \ c I ^•''^CO [m^s] 

(31) 

The Planck spectral response conversion coefficients are 
therefore given by the following ratios 



^KCMB_to_MJy/s] 



^KCMB_to_ySZ = 



dv T(v)b[, 



X 10 



,20 



dv t(v)(vc/v) 



MJy/sr 



KcMB 



(32) 



dv T(v)b[, 



" kT I exp [hv/(kT)\ - 1 



CMB 



[KcMB 



(33) 



^ The variable v is used to denote velocity in units of km s ^ , while 
the variable y is used to denote frequency in units of Hz (or equivalent). 



Jdv t(v)(vc/v) 
^dvT{y)(ylv,r 

Jdv t(v)(vc/v) 



Hz 
Hz 



dVT(v)(v/Vcf 



and 



Uco = 



rtvco)(^2°)i.; 



Byiy, T) 



v,^By(vc,T) 



Hz 
Hz 



(34) 



(35) 



RJI 



ro 



dVT(v)K 



KcMB 



Krj km s 



(36) 



3.2.3. HFI Unit Conversion and Colour Correction 
Coefficients 

This section presents unit conversion and colour correction co- 
efficients resultant from the above relations and HFI detector 
spectra. Similar values for the LFI may be found in Planck 
Collaboration II (2013, e.g.. Table 8). In addition to processing 
data from individual detectors, HFI data are processed to pro- 
vide band-average frequency maps, and sub-band average fre- 
quency maps. The sub-band frequency maps are comprised of 
three sets, DetSetl, DetSet2, and SWB (i.e., spiderweb bolome- 
ters only). These groupings for the HFI detectors are summa- 
rized in Table 4. Results presented here are restricted to the av- 
erage spectra. Similar data for individual detectors is available in 
Planck Collaboration ES (2013). Table 5 provides multiplicative 
unit conversion coefficients to go from MJy sr"^ to K5 brightness 
temperature. ^ Table 6 provides unit conversion coefficients for 
band-average spectra. An example dust colour correction coeffi- 
cient is also provided in the table, which has assumed a dust SED 
with a spectral index of 4. Figure 13 illustrates the variation in a 
colour correction coefficient as a function of power-law spectral 
index. Colour correction coefficients for planets within our Solar 
system (HFI calibration sources) have also been determined; 
Table 7 provides these coefficients for the band-average spec- 
tra; these are needed, e.g., for the HFI beam calibration (Planck 
Collaboration VII 2013). The CO conversion coefficients are 
provided in Sect. 4.2.1 and in Planck Collaboration ES (2013). 

The multiplicative unit conversion and colour correction co- 
efficients are to be used as follows. Take, for example, a dust 
region within a Planck 100 GHz band-average map with an es- 
timated intensity of IOKcmb- To convert this intensity to an 
equivalent specific intensity in MJy sr"^ (IRAS), the original in- 
tensity should be multiplied by the unit conversion coefficient 
of 244.1 (see Table 6) to obtain a brightness intensity of 2441 
MJy sr"^ If the dust can be approximated to follow a SED spec- 
tral profile with a spectral index of 4, then the colour correc- 
tion would be applied by multiplying the 2441 MJy sr"^ by the 
colour correction coefficient of 0.8938 to convert /r(9m a spec- 
tral index of -1 to 3. spectral index of 4; yielding 2182 MJy sr"^ 
The fourth column of Table 6 provides the combined unit con- 
version / colour correction as the product of the second and third 
columns. If the same IOKcmb intensity were found within a 
100 GHz detsetl map, then the resulting colour-corrected dust 



^ This unit conversion is independent of spectral index, and the 
brightness temperature does not imply a Rayleigh- Jeans spectral pro- 
file. 
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Table 4. HFI sub-band average groupings 



Band [GHz] 


DetSetl 


DetSet2 


SWB 


100 GHz 


la/b, 4a/b 


2a/b, 3a/b 




143 GHz 


la/b, 3a/b 


2a/b, 4a/b 


5,6,7 


217 GHz 


5a/b, 7a/b 


6a/b, 8a/b 


1,2,3,4 


353 GHz 


3a/b, 5a/b 


4a/b, 6a/b 


1,2,6,7 


545 GHz 


1,2 


4 




857 GHz 


1,2 


3,4 






-6 -4 -2 2 4 6 

Powerlaw Spectral Index 



Fig. 13. Colour correction coefficients for varying spectral index 
for the HFI band-average spectra. Similar data for individual de- 
tectors and DetSet subsets are available in Planck Collaboration 
ES (2013). 

Table 5. HFI flux density to brightness temperature unit conver- 
sion coefficients 



Band [GHz] 


Uc [KRjMJy-i sr] 


100 


0.0032548074 


143 


0.0015916707 


217 


0.00069120334 


353 


0.00026120163 


545 


0.00010958025 


857 


0.000044316316 



intensity would be 2176 MJy sr"^ . To colour correct from a spec- 
tral index of 4 to 3. spectral index of -1, one would divide the 
intensity by the colour correction coefficient found in the third 
column of Table 6. 

To demonstrate the relative contribution of the various re- 
gions of the spectral bands on the CMB dipole signal, Fig. 14 
illustrates the b[,T(y) product for the band-average spectra (sim- 
ilar results are found for the individual detector spectra). It is 
important to note that for the 857 GHz channel, and partly for 
the 545 GHz channel, the unit conversion integral (Eq. 8 and 
32) is dominated by the low-frequency portion of the band. Any 
residual systematics in the transmission spectra may cause un- 
desired errors, e.g., in the conversion of 857 GHz (or 545 GHz) 
data to units of Kcmb • H is preferred, for example, to convert the 
lower frequency channel data to MJy sr"^ when needed, rather 
than converting the 857 GHz data into Kcmb units. For compar- 



Wavenumber [cm ""l 
5 10 20 30 




Frequency [GHz] 



Fig. 14. Product of the the HFI band- average spectra, T(y), with 
the CMB dipole spectral profile, i.e., /?^T(y) (Eq. 8) shown using 
solid curves. The nominal spectra shown for as dotted curves for 
reference. 



Wavenumber [cm ""l 




Frequency [GHz] 

Fig. 15. Product of the band-average spectra with an example 
dust spectrum (T^ =18K, =1.5, Eq. 27) shown using solid 
curves. The nominal spectra are shown as dotted curves for ref- 
erence. 



ison, Fig. 15 illustrates the product of the band-average spectra 
with a sample dust spectrum, with a modified blackbody of dust 
temperature, =18 K andySj =1.5 (see Eq. 27). 

Another investigation was conducted in order to understand 
the eff'ect of a hypothetical systematic bias in the spectral re- 
sponse data. This study involved scaling the transmission spec- 
tra by an scaling term with a linear dependence on frequency, 
such that it is unity valued on the nominal band centre (i.e., 100, 
143, 217, 353, 545, and 857 GHz for the respective bands), is 
centred on the nominal band centre frequency, and has a linear 
deviation towards a specified value, mj-ef, at the Vc ± 15 % band 
edges. Fig. 16 illustrates the normalized variation of a combined 
unit conversion and colour correction for each of the HFI bands, 
over a range of linear slopes spanning m^ef ^ [-2 %, 2 %]. The 
selected example illustrates a conversion from Kcmb to MJy sr"^ 
and a colour correction from spectral index of -1 to the dust 
profile described above. The spectral uniformity of the reference 
bolometer used in characterizing the HFI detector spectral re- 
sponse (see Planck Collaboration ES (2013)) is estimated to be at 
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Table 6. HFI unit conversion coefficients 



Band 


Uq IRAS 


dust 


(^C TRAS)(^r dust) 


^CSZ 


[GHz] 


[MJy/sr/KcMB] 






[MJy/sr/KcMB] 


[Ysz/Kcmb] 


1 VJVJdVg 


244.1 +0.3 


0.8938 


+ 0.0019 


918 9 + 


-0.24815 


+ 0.00007 


1 OO-Hpt^iPtl 

1 VJVJ-LlCldCLl 


244.9 + 0.4 


0.889 


+ 0.003 


0\1 fs +04 


-0.24833 


+ 0.00010 






0.896 


+ 0.003 


218.4 + 0.4 


-0.24807 


+ 0.00009 


143-avg 


371.74 + 0.07 


0.9632 


+ 0.0004 


358.04 + 0.07 


-0.35923 


+ 0.00006 


1 zL^-HptQPtl 


^^SS 0^ + 1 S 


1.0058 


+ 0.0009 


i,fj-l 1 < + 1 S 


-0.35398 


+ 0.00011 


1 4^-rlpt<iPt9 

1 T- J Lie 1 5> CIZ/ 




0.9773 


+ 0.0008 


Q(^n +01^ 


-0.35743 


+ 0.00010 




'57c <Q +0 14 

J/O.JO X VJ.IH- 


0.9238 


+ 0.0008 


^4Q 74 +014 


-0.36446 


+ 0.00011 


217-avg 


/loo /' c\r\ r\ r\ \ ^ 

483.690 + 0.012 


0.85895 


+ 0.00011 


415.465 + 0.012 


5.152 


+ 0.006 


9 1 7-rlptQPt1 


480 + 09 


0.88411 


+ 0.00016 


AT) A f^Q + 09 


7.212 


+ 0.019 


917 rlptcpt9 


480 ^1 4 + 01 Q 


0.88235 


+ 0.00017 


49^ 804 + 01 Q 


7.046 


+ 0.018 


217-SWBs 


486.331 +0.018 




+ U.UUUiJ 


408.855 + 0.018 


4.236 


+ 0.006 


353-avg 


987 4S0 + OOQ 


0.85769 


+ 0.00011 


94fS ^4^ + OOQ 


0.161098 


+ 0.000011 


^^^-HptQPtl 

J J J LlCloCll 


98Q 690 + 01 9 


U.ooZj J 


+ U.UUUiJ 


9SS 606 + 01 9 


0.163757 


+ 0.000014 


353-detset2 


287.967 + 0.013 


0.86548 


+ 0.00014 


249.229 + 0.013 


0.160904 


+ 0.000014 


353-SWBs 


286.786 + 0.011 


0.84997 


+ 0.00014 


243.759 + 0.011 


0.160456 


+ 0.000013 


545-avg 


58.04 +0.03 


0.85444 


+ 0.00016 


49.59 +0.03 


0.06918 


+ 0.00003 


545-detsetl 


58.02 +0.03 


0.8513 


+ 0.0002 


49.39 +0.03 


0.06924 


+ 0.00004 


545-detset2 


58.06 +0.05 


0.8612 


+ 0.0003 


50.00 +0.05 


0.06905 


+ 0.00005 


857-avg 


2.27 +0.03 


0.9276 


+ 0.0002 


2.09 +0.03 


0.0380 


+ 0.0004 


857-detsetl 


2.26 +0.03 


0.9231 


+ 0.0003 


2.08 +0.03 


0.0380 


+ 0.0005 


857-detset2 


2.27 +0.04 


0.9333 


+ 0.0003 


2.11 +0.04 


0.0380 


+ 0.0005 




°-2 -10 12 



Slope [%] 

Fig. 16. Variation of combined unit conversion / colour correc- 
tion with a hypothetical systematic spectral error as described in 
the text. 



the level of 1 %; this is the motivation behind the type of system- 
atic distortion introduced in this study. The figure demonstrates 
that a systematic reference spectral flatness error, as described 
above, results in biased unit conversion / colour correction co- 
efficients. The introduced coefficient bias has a magnitude that 
scales linearly with the slope of the spectral flatness systematic 
error introduced. The 857 GHz channel fluctuates the most due 
to the the dominance of the unit conversion on the low-frequency 
region of the band, and the dominance of the selected colour cor- 
rection on the high-frequency portion of the band, as illustrated 
in Figs. 14 and 15. 



3.2.4. Unit Conversion and Colour Correction Software Tools 

A set of software tool, written in the Interactive Data Language 
(IDL) has been developed for distribution with the HFI detec- 
tor spectra and Planck data. This tool package, herein the UcCC 
package, uses the transmission spectra provided in a .fits file 
format, and computes the unit conversion and colour correction 
factors using the relations derived above. The UcCC tools may be 
used to determine colour corrections for a variety of spectral pro- 
files, including powerlaw and modified blackbody as described 
above as well as user defined source spectra, these tools also pro- 
vide coefficient uncertainty as an optional output. The UcCC code 
package may be obtained from the Planck Legacy archive, with 
further details on their use provided in Planck Collaboration ES 
(2013). 



3.3. Error Propagation 

The uncertainty of the HFI detector spectral response and band- 
average spectral response was propagated to the coefficient fac- 
tors described above. For each of the correction factors shown 
(Eq. 32 - 36), the measured spectral transmission profile, T(y), 
is used along with the associated spectral uncertainty. An addi- 
tional uncorrelated uncertainty was repeatedly introduced to the 
spectrum under consideration, where a Gaussian noise distribu- 
tion in frequency space was weighted by the uncertainty in T(y), 
with the resultant spectrum plus noise then used in the calcula- 
tion of the desired coefficient. The uncertainty values for coef- 
ficients shown in this work correspond to the statistical evalua- 
tion of 10000 trials in each instance. These uncertainties reflect 
only the propagation of the spectral uncertainty into the coeffi- 
cient and may not represent the full uncertainty in every case. An 
evaluation of some potential sources of systematic uncertainty is 
included in Sect. 4 below. 
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Table 7. HFI planet colour correction coefficients 



Band [GHz] 


Mars 


Jupiter 


Saturn 


Uranus 


Neptune 


100-avg 
143-avg 
217-avg 
353-avg 
545-avg 
857-avg 


0.9613 +0.0017 
1.0069 +0.0009 
0.9355 +0.0002 
0.93365 + 0.00005 
0.93603 + 0.00010 
0.98273 + 0.00014 


0.962 + 0.005 
1.007 +0.004 
0.9376 +0.0016 
0.9385 +0.0004 
0.92946 + 0.00009 
0.99918 + 0.00014 


0.963 + 0.004 
1.008 +0.003 
0.9529 +0.0010 
0.94145 + 0.00014 
0.81721 + 0.00009 
0.99923 + 0.00014 


0.9692 +0.0011 
1.0122 +0.0002 
0.94897 + 0.00004 
0.94870 + 0.00004 
0.95148 + 0.00009 
0.99210 + 0.00013 


0.9741 +0.0010 
1.0127 +0.0002 
0.96195 + 0.00004 
0.93676 + 0.00004 
0.96409 + 0.00009 
0.99815 + 0.00013 



3.4. Band-average Coefficient Maps 

Using the Wmj(6, (p) maps described above (Eq. 6), band-average 
spectra may be computed for every pixel of a sky map (rather 
than integrating the relative weights across the map). Thus, for 
frequency channel maps, coefficient sky maps may be generated 
for individual surveys and combinations of surveys. For individ- 
ual detectors, the response is expected to be constant across the 
sky, but for channel average data, the sky coverage and relative 
noise causes variations in the proportional averaging. An exam- 
ple band-average coefficient map is shown in Fig. 17, where the 
353 GHz combined Kcmb to MJysr"^ (IRAS) unit conversion 
and IRAS to dust (see Sect. 4.2.2) colour corrections are shown 
for the nominal and individual surveys. 

These coefficient maps can be used to investigate differences 
between surveys, and the coefficient variation across a region of 
interest can be compared to the magnitude of other sources of 
error, providing a probe of the effect of scan strategy, integra- 
tion time, and relative intra-band detector noise levels on map 
consistency and unit conversion coefficients. Histograms of the 
coefficient sky maps provide a verification of the Monte-carlo 
uncertainty estimates in much the same way that the weight fac- 
tor histograms verify the band-average scaling factors (see Fig. 
12). For coefficient map histogram distributions that are rela- 
tively narrow, i.e., of the same width as the Monte-Carlo simu- 
lation based uncertainty estimates, data processing and analysis 
may be simplified by using static conversion coefficients as op- 
posed to coefficient maps. While these maps are not being pro- 
vided under the current data release, the discussed histograms 
have been used to validate the use of scalar unit conversion and 
colour correction coefficients as opposed to the discussed coef- 
ficient maps. The maps within Fig. 17 demonstrate the variation 
of the band-average unit conversion / colour correction coeffi- 
cient to be small, typically less than 0.1 %. Histograms of these 
maps are shown below in Fig. 27. This discussion will continue 
in greater detail within future work as these levels of coefficient 
fluctuations at the map level are expected to become more im- 
portant with polarization analysis. 




Fig. 17. 353 GHz band- average dust unit conversion and colour 
correction coefficient maps. The maps show the deviation of 
the coefficient about the median. The medians in each case are 
278.970 (a.), 278.960 (b.), and 278.970 (c.) MJy sr'^ K"!^ 



4. Discussion 

This section evaluates the HFI detector and band-average spec- 
tral response data, and the associated unit conversion and colour 
correction algorithms and coefficients. This is done through 
comparison of these data with HFI flight data. HFI detector 
bandpass mismatch, i.e., the relative difference between an in- 
dividual detector spectrum and its band-average counterpart, 
can be compared with variations between individual detector 
and channel average detector results. Examples include CO, SZ, 
dust, etc. Out of band signal rejection is verified through com- 
parisons with the HFI Zodiacal light observations, where any 



out-of-band sky signal received would be evident in the data. 
As the progress of the HFI polarization data advances, compar- 
isons between spectral response predicted polarization leakage 
and observed intensity to polarization leakage in HFI maps can 
also be used to verify the accuracy of the HFI spectral response 
data. 
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4.1. Out-of-band Signal Rejection 

Each of the HFI bands has a filter stack of 5 low-pass filters 
with varying cut-off' frequencies in order to achieve suitable out- 
of-band signal rejection. In addition, the 545 and 857 bands 
also have a high-pass filter dictating the cut-on frequency, which 
also serves to provide some rejection at much higher frequen- 
cies. When the transmission of the individual filters comprising 
each HFI filter stack was measured, a series of FTS scans were 
recorded with external low pass filters up to 650 cm"^ (approx. 
19.5 THz, 15.4 Jim). This was not done as a single measurement 
but a series of experiments to allow the transmission in-band to 
be measured to higher resolution and S/N. For each of the low 
pass-filters, the transmission at the high-frequency end dropped 
to about 10"^ to 10"^ beyond the cut-off' frequency. Given that 
there are 5 filters in a stack, the out-of-band rejection is very 
high. The spectral transmission for all frequencies greater than 
3 THz is less than 10"^^ for all HFI bands, reaching as low as 
10"^^ in some cases. This level of out-of-band signal 
rejection is confirmed by the Zodiacal light data described be- 
low. 

4.1 .1 . Zodiacal Light Verification 

As mentioned above, optical filters are used to prevent out-of- 
band light from impinging on the HFI detectors and registering 
as signal. While these filters are thoroughly characterized and 
tested before ffight, it is of interest to try to confirm, as well as 
possible, their behavior in situ. In this section, the HFI observa- 
tions of the Zodiacal light, specifically 100 GHz data, are used 
to place an upper limit on out-of-band spectral contribution, i.e., 
spectral leaks. The spectral profile of Zodiacal emission is well 
understood (e.g., Kelsall et al. (1998)), with Zodiacal emission 
much brighter in the mid-infrared than in the millimeter wave- 
length range and very little Zodiacal-correlated signal observed 
at 100 GHz (Fixsen & Dwek 2002). While the Diff'use Zodiacal 
Cloud is observed in the higher- frequency HFI channels (Planck 
Collaboration XIV 2013), it is not detected within the lower fre- 
quency HFI channels. A model Zodiacal light spectral emission 
profile is shown in Figure 18, with a dashed line representing 
the reduced emissivity for frequencies less than 2 000 GHz, as 
expected. 

To set limits on a hypothetical leak at short wavelengths, 
the curve in figure 1 8 is frequency integrated to find the ex- 
pected Zodiacal cloud emission ffux density. To estimate the 
expected in-band ffux density, an in-band emissivity of approx- 
imately 0.041 (see Planck Coffaboration XIV (2013)) is cou- 
pled with the Zodiacal emission proffle at an eff'ective black- 
body temperature of 286 K (Zioo ghz below). This is combined 
with the nominal 100 GHz channel bandwidth (Avioo ghz below, 
approximately 33 GHz), where the product yields the expected 
in-band Zodiacal signal. This product is then divided by the to- 
tal Zodiacal ffux density estimated by the frequency integration 
of the Zodiacal spectral proffle. Thus, an upper limit on any 100 
GHz high frequency spectral leakage, Lu, is given by: 

^ ^100 GHz • Avioo GHz 

X ^-^W . (37) 
< 1.2-10"^ 

The frequency range used in the above integration is 3 THz - 
150 THz (equivalent to lOOyum - 2yum). This integration region 
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Fig. 18. The spectrum of the diff'use cloud in the Ecliptic plane. 
The blue curve represents the spectrum one would see assuming 
that uniform cloud emissivity of unity. The dotted red line shows 
the modiff cation expected with a l/A^ emissivity proportionality 
for frequencies below 2 000 GHz (i.e., wavelengths longer than 
150yum). 



covers the dominant signal and is well above the 100 GHz band 
cut-off'. Diff'erences in the integration between this range and 
larger ones result in changes to the spectral leak upper limit at 
the sub-percent level, and thus do not eff'ect the calculation at the 
desired accuracy. 

The lack of any Planck detection of the Diff'use Zodiacal 
cloud signal at 100 GHz (Planck Collaboration XIV (2013)), in- 
dicating high out-of-band Zodiacal signal rejection, thus places 
conffdence in the out-of-band signal discrimination of the HFI 
detector optical fflters. Through this analysis, out-of-band dis- 
crimination has been shown to be greater than 10^ for the HFI 
100 GHz detectors. Less stringent limits can be set on smaller 
subsets of the band including wavelengths greater than 100/im. 

4.2. Bandpass Missmatch 

Variations of the spectral response of individual detectors from 
the band-average response can be compared with similar varia- 
tions in HFI data. An example of this is the observation of CO 
emission within some of the HFI bands. As CO emission is in- 
trinsically narrow-banded, diff'erences in CO sensitivity for dif- 
ferent HFI detectors are easily compared to the HFI spectral re- 
sponse variations. Other components of the sky signal may also 
be used to perform similar comparisons (e.g., SZ, polarization 
leakage, dust, etc.). 

Polarization leakage is an eff'ect where spectral mismatch be- 
tween detectors in the same frequency band induces a reduc- 
tion/loss in an unpolarized intensity map and an over-estimation 
of the corresponding polarization map. This eff'ect, which is a 
result of the assumptions required within the mapmaking algo- 
rithms, can also cause polarization to intensity leakage, but this 
is of less concern, however, due to relatively weak amplitudes of 
polarization signals coupled with the low levels of leakage ex- 
pected. While a discussion of the polarization aspects of the HFI 
spectral response will be withheld until the release of Planck po- 
larization data, the concept is introduced here to provide context 
to the current work as the same eff'ect that causes polarization 
leakage allows a sky-based estimate of the respective colour cor- 
rection / unit conversion coefficients (see Sect. 4.2.2). 
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4.2.1 . CO Bandpass Verification 

Analysis of HFI data, using component separation methods that 
include a CO emission component (Planck Collaboration XIII 
(2013) and Planck Collaboration XII (2013)), has provided all- 
sky maps of CO emission for the first three CO rotational tran- 
sitions, i.e., CO /=1^0, CO /=2^1, and CO /=3^2. Select 
molecular cloud regions with well-known CO emission prop- 
erties provide external validation of the HFI CO maps. One of 
the validation observations used was the DAME survey (Dame 
et al. 2001), which observed the CO /=1^0 transition. As the 
HFI CO maps are the result of component separation performed 
on maps in Kcmb units, and the external CO observations are 
typically available in units of velocity integrated brightness tem- 
perature, i.e., Kkms"\ this allows for a sky-based estimate of 
the CO unit conversion coefficients from units of Kkms"^ to 
units of Kcmb • These coefficient estimates are then compared 
against the bandpass-based coeflftcients, i.e., those based on the 
pre-flight measured spectral transmission data and the respective 
unit conversion relation above (Eq. 36). Another benefit of com- 
paring the bandpass and sky -based CO results is the potential for 
improved understanding of systematic uncertainties in the band- 
pass data. 

While a detailed discussion regarding the derivation of the 
sky-based CO unit conversion coefficients is found elsewhere 
(Planck Collaboration XIII 2013), this work discusses only those 
details relevant to the comparison of the sky and bandpass CO 
coefficients. The bandpass CO coefficients provide a direct con- 
version from diff'erential CMB emission to that of CO, with the 
caveat that a precise knowledge of the CO unit conversion co- 
efficients from the bandpass data is difficult to obtain in that 
the spectral resolution of the bandpass data is much broader 
than any observed CO emission features. Thus, the bandpass 
CO coefl&cients are an estimate based upon under-resolved spec- 
troscopic measurements. The sky-based CO coefficients, on the 
other hand, provide a relative conversion coefficient based on 
the variation in spectral transmission between detectors within a 
common frequency band. 

The linear combination method of CO signal extraction used 
to obtain the HFI CO products involves using a weighted sum 
derived to maximize the contrast between the desired component 
and its residuals (see Hurler et al. (2010) for details). Within the 
use of this method in the CO extraction, the DAME data set is 
used as a calibration template for the weighted map sums within 
a frequency band (i.e., the Planck 100, 217, or 353 GHz bands). 

As demonstrated in Eqs. 1 and 5 of Planck Collaboration 
XIII (2013), the weighted coefficient for an individual detector 
signal is given by w/F/, where w/ is the MILCA relative weight, 
and Fi is the CO unit conversion coefficient. For the sky-based 
CO coefficients, the coefficient estimate is based upon the corre- 
lation of the CO component separation output with the DAME 
survey; thus, the relative CO and CMB weighting of the band- 
pass spectra (i.e., the numerator and denominator of Eq. 36) are 
determined indirectly, without the bandpass data. To illustrate 
the diff'erence between the relative transmission at the CO rota- 
tional transition frequencies (e.g., just the numerator of Eq. 36), 
and the relative unit conversion for the CO transitions (e.g., all of 
Eq. 36), Fig. 19 compares these two sets of parameters. Although 
the two relative values demonstrate a correlation, they do not 
demonstrate perfect agreement. It is therefore important that CO 
unit conversion coeflftcients account for both the relative diff'er- 
ences in spectral transmission at the CO rotational transition fre- 
quencies as well as the relative spectral transmission variations 
within the CMB signal over the entire band of a given detector. 
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Fig. 19. A comparison of the relative spectral transmission at 
the CO transition frequency with the relative CO unit conversion 
coefficient for the 100 (a.), 217 (b.), and 353 (c.) GHz bands. 

Through modeling of CO emission in a molecular cloud as 
a function of CO rotation temperature, it can be shown that the 
relative line intensities vary with CO rotation temperature, where 
the lower order transitions dominate for very low temperatures 
(e.g., less than 10 K). As the relative line intensity for, e.g., the 
CO /=2^1 and CO /=1^0 transitions varies with tempera- 
ture, and the lowest CO transition is used in the HFI/DAME CO 
map calibration, the uncertainty of the HFI CO maps is expected 
to vary with CO temperature. This is consistent with the uncer- 
tainty varying with the observed CO intensity. Thus, in compar- 
ing the sky and bandpass CO coefficients, it is important to note 
that the relative intensities of CO emission change with CO ro- 
tational temperature and that for the higher order rotational tran- 
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Fig. 20. Comparison of the relative CO specific intensity for the 
first three rotational transitions, for a varied CO rotation temper- 
ature. The normalization employed removes any dependence on 
the column density. 



sitions, the sky CO coefficients are a generalized scalar value 
representing a dynamic quantity. 

The sky-based coefficients do not provide a conversion that 
includes an absolute calibration, but represent the relative lev- 
els between detectors of the same frequency band; the absolute 
calibration comes from the external data sources used for valida- 
tion. This is a result of the linear combinations used in combin- 
ing detector signals in such a way as to enhance certain spectral 
components while reducing other spectral contaminants (Hurler 
et al. 2010). 

Another diff'erence between the sky/DAME and bandpass 
CO coefficients is the isotopic content of the calibration. The 
sky -based coefficients provide an estimate of CO emission from 
all isotopic species, while the bandpass CO coefficients are cal- 
culated for isotopic contributions individually (e.g., the emis- 
sion frequencies vary with isotope). The bandpass coefficients 
for individual isotopes could be combined provided information 
regarding the isotopic ratios were available for a given region. 
Thus, there is a normalization/rescaling step needed to compare 
the two values, and a complete comparison is unconstrained at 
present without additional calibration data. 

Additionally, the sky-based coefficients are based on unpo- 
larized intensity, so the PSB detector pairs are combined in the 
form {A + B)I2 io reduce the influence of any polarized signal. 
Any uncertainty in the data introduced as a result of this step 
must be propagated through to the sky CO coefficients. 

For these reasons the sky and bandpass CO coefficients are 
expected to show a correlation, but are not expected to exhibit 
perfect agreement. This is particularly true of an absolute com- 
parison; relative diff'erences should be proportionate, but the 
absolute values of the bandpass and sky coefficients are not 
expected to be directly compatible (Planck Collaboration XIII 
2013). 

Eff'orts were taken to understand the accuracy of the band- 
pass CO coefficients in light of the diff'erences observed between 
the sky and bandpass CO coefficients. As demonstrated in Eq. 
36, the relative spectral transmission over a very narrow fre- 
quency range, in principle a delta function, is needed in order 
to obtain a CO unit conversion coefficient. This ideally requires 
knowledge of the HFI spectral response to much finer spectral 
resolution than is available. Each of the other unit conversion and 
colour correction calculations are based on relatively broad spec- 



tral features that are not dominated by uncertainties within just 
one spectral bin. Thus, it is difficult to ensure that the nominal 
transmission and uncertainty for a single spectral bin are consis- 
tent with the values provided within the detector bandpass data. 
In order to investigate the eff'ects of CO transmission interpola- 
tion and uncertainty errors on the CO conversion coefficients, the 
bandpass CO coefficients were repeatedly calculated after vary- 
ing levels of smoothing were applied to the data. Smoothing was 
done in steps of one from one to nine. The smoothing results in 
noise averaging within an increasingly broad spectral bin width, 
so the intrinsic spectral transmission uncertainty is reduced, but 
the interpolation error may be increased. Fig. 21 illustrates the 
changes in CO coefficients upon introduction of this spectral 
smoothing. It is evident that the spread in coefficient values with 
smoothing factor is greater than the coefficient uncertainty in 
many cases. Thus, the bandpass CO coefficients have been re- 
vised to include this result. Within the Planck Collaboration XIII 
(2013) paper, this is reflected as an increase in the error bars of 
the bandpass CO coefl&cients. 

Another check performed on the bandpass CO coefl&cients 
was the introduction of a linear slope to the detector spectra, as 
described in Sect. 3.2.3 and demonstrated in Fig. 16 for a com- 
bined unit conversion and dust colour correction. As was found 
to be the case above, the addition of a linear slope to the spectra 
resulted in a linear change in the CO coefficients. In this instance 
the changes observed in the CO coefficients were well within the 
quoted uncertainty, so the uncertainty is not underestimated in 
this respect. Fig. 22 illustrates the shift in CO coefficients caused 
by the introduction of a linear scaling of the spectra. 

The variation between the sky-based and bandpass-based CO 
unit conversion coefficients is demonstrated in Figs. 23 - 24 for 
the HFI 100, 217, and 353 GHz channels. The first figure in- 
dicates the correlation between the normalized sky and band- 
pass coefficients, while the second shows the sky and bandpass 
coefficients grouped together about their respective feed horns. 
Plotting the sky coefficients on this scale required scaling the 
relative sky coefficients about the bandpass coeflftcient mean. 
This was done to allow comparison with the ^^CO coefl&cients 
as well as the ^^CO coefficients (Fig. 24 also includes the ^^CO 
coeflftcients although these were not used in the sky coeflftcient 
scaling). The ^^CO coefficients are expected to have a stronger 
influence on the 100 GHz channel than for the higher frequency 
channels (Planck Collaboration XIII 2013). The plots indicate 
the detector type as either the individual PSB detectors (A and 
B separately), the combined PSB detectors ((A -h B)/2), or the 
spiderweb (SWB) detectors. Although this comparison is under- 
constrained, and the two sets of coefficients do not represent ex- 
actly the same quantity, as has been discussed above, there is a 
general agreement between the sky-based and bandpass-based 
CO coefficients. The CO analysis also required a dust template 
as part of the component separation. Thus, a sky-based estimate 
of HFI dust colour correction coefficients was a by-product of 
this analysis. There is excellent agreement found between the 
bandpass- and CO/sky -based dust colour correction coefficients; 
more details on this are provided in Sect. 4.2.2, with the results 
presented in Fig. 28. 

There remain diff'erences in the sky- and bandpass-based CO 
unit conversion coeflftcients. The analysis of CO data in the 2013 
Planck data release is based upon the sky-based CO coefficients. 
The diff'erences between the two approaches are acceptable for 
present analysis, but are increasingly important for future work, 
including the analysis of Planck polarization data. Thus, while 
there is a correlation between the two approaches, a better under- 
standing of these diff'erences is required to gain a deeper under- 
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Fig. 21. Comparison of the CO unit conversion coefficients with 
a smoothing / re-binning factor applied to the spectra going in to 
the unit conversion algorithm, for the 100 (a.), 217 (b.), and 353 
(c.) GHz detectors. 



Fig. 22. Comparison of CO unit conversion coefficients with an 
additional linear scaling of the input spectra, for the 100 (a.), 217 
(b.), and 353 (c.) GHz detectors. 



standing of the Planck data. This comparison presents the cur- 
rent standing of a work in progress. 

4.2.2. Dust Colour Correction Bandpass Verification 

As a result of work investigating polarization leakage, a study 
was conducted to estimate the on-sky integration of dust SEDs 
for dusty regions directly using flight data, i.e., without using 
the spectral transmission data. This was done to investigate the 
compatibility of the two methods. 



To derive the dust spectral mismatch for a given HFI fre- 
quency band, dust colour correction coefficients for each individ- 
ual detector must be known for the bolometers/detectors within 
the given band. Difl'erences between the dust content of individ- 
ual bolometer maps and the corresponding frequency channel 
map can be used to estimate these coefficients. In performing 
this analysis, care must be taken to understand all of the spec- 
tral components included in the maps (e.g., CMB, free-free, CO, 
etc.), and calibration errors, beam errors and beam difl'erences, 
and any residual polarization signal within the bolometer inten- 
sity maps. 
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Fig. 23. Relationship between the CO unit conversion coeffi- 
cients based on the CO- sky maps and the spectral response data 
for the CO /=1^0 transition within the 100 GHz HFI band 
(top), the CO J=2^\ transition within the 217 GHz HFI band 
(middle), and the CO /=3^2 transition within the 353 GHz HFI 
band (bottom). 
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Fig. 24. Comparison of CO unit conversion coefficients, scaled 
about the bandpass coefficient amplitudes for the sky and band- 
pass coefficients. This scaling allows comparison with both the 
^^CO and ^^CO bandpass coefficients; the bandpass isotopic co- 
efficients are horizontally offset in opposing directions for clar- 
ity. 



To avoid polarization bias, the study was restricted to regions 
surrounding the ecliptic poles, where there is a large variety in 
the crossing angles for multi-scan observations of a given re- 
gion with a single detector. Stokes /, 2, and U maps were pro- 
duced for individual detectors in addition to the standard /, 2, 
and U frequency channel maps nominally produced by the data 
pipelines. 

From the ecliptic polar maps, there were three regions se- 
lected for this study, all specifically selected for their relative 



component/content diff'erences. The Large Magellanic Cloud 
(LMC, (/,/?)=(280?5, -32?8), Planck Collaboration XVII (2011)) 
was selected as a familiar target with known dust properties. 
Another source selected was a region slightly North of the LMC 
in the Chamealeon constellation (CHA, (/,/?)=(287?6, -24? 1)), 
and the final region used in this study was near the North Ecliptic 
Pole (NEP, (/,/?)=(103?0, 18?8)). For each of the HFI frequency 
bands, and all of the individual detectors, square maps of 13° 
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Fig. 25. 353 GHz maps of the ecliptic polar maps for the LMC 
(left), CHA (middle), and NEP (right) regions of the sky. 



width, with 4' pixel resolution, were extracted from the /, Q, 
and U maps described above. The dust SED properties of the 
LMC from a previous Planck publication are = 21.0 ± 1.9 K 
and = 1.48 ± 0.25 (Planck Collaboration XVII (2011)); as 
the sky -based dust colour correction coefficients for all sources 
remain consistent (see Fig. 28), the LMC dust properties and 
uncertainties mentioned above are used to determine the band- 
pass coefficients included in this section of this work, along with 
their respective uncertainty. Fig. 25 provides the HFI observa- 
tions over these regions at 353 GHz. The corresponding band- 
average unit conversion / colour correction coefficient maps (see 
Sect. 3.4) for each band were used to verify the uniformity of 
the expected band-average colour correction across these regions 
of the sky. While Fig. 17 provides examples of full-sky band- 
average coefficient maps. Fig. 26 provides a similar example 
over the three specified sky regions used in the sky/bandpass dust 
colour correction coefficient comparison. The changing level of 
coverage between surveys is especially apparent for the NEP re- 
gion. Histograms of these band-average coefficients for the full 
sky and dusty regions are provided in Fig. 27. 

A component separation is performed to isolate CMB, free- 
free, CO, and dust emission from the 100, 143, 217, and 353 
GHz frequency channels. Regions of high free-free and high CO 
emission are masked in the dust maps, leaving residual maps 
primarily containing CMB and dust. The CMB signal is re- 
moved by subtracting the frequency-band map from the individ- 
ual detector maps as described below. Let U d represents the dust 
colour correction for the frequency-band map of a given chan- 
nel, and let U d i represent the individual bolometer dust colour 
correction coefficient. The intensity map, /, for the channel will 
thus be given by / = C -I- UdD, where C is the CMB compo- 
nent, and D is the dust component (note: C and D represent the 
actual observable intensity, not just what the detector sees). The 
detector intensity map will be given by // = C -\- Ud iD. The 
regression coefficient resultant from the correlation of the diff'er- 
ence map // - / with the channel map / can be used to derive the 
relative dust colour correction coefficient as follows: 
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Thus the relative dust colour correction coefficient may be deter- 
mined for each individual HFI detector by employing the above 
relation. It is important to note that the dust colour correction 
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Fig. 26. 353 GHz maps of dust band- average unit conversion / 
colour correction (see Sect. 3.4). the maps show the deviation in 
the combined unit conversion / colour correction from the full- 
sky map median values. The full- sky median coefficient values, 
i.e., off'sets, used for the Full (a.,b.,c.). Nominal (d.,e.,f.). First 
(g.,h.,i.), and Second G.,k.,l.) surveys are, 278.980, 278.970, 
278.960, and 278.970 MJy sr'^ Kcmb, respectively. The figure 
columns represent the ecliptic polar maps for the LMC (a., d., g., 
and j.), CHA (b., e., h., and k.), and NEP (c, f., i., and 1.) regions 
of the sky. 



coefficients in the above expression convert from the dust spec- 
tral profile to the CMB spectral profile; in the notation of the 
HFI unit conversion and colour correction syntax, this is equiv- 
alent to the inverse (i.e., U^) of a colour correction from the 
dust spectral profile to a powerlaw spectral index of - 1 followed 
by a unit conversion from a powerlaw spectral index of -1 to 
diff'erential CMB temperature units. 

Fig. 28 illustrates the excellent agreement between the sky- 
based Dust colour correction coefficients and those based on the 
HFI spectral response (i.e., Eq. 32 and 35) for the 100, 217, and 
353 GHz spectral bands. This level of agreement is not found 
amongst the 143 GHz detectors, however. The source of this dis- 
crepancy is under investigation. It should be noted that the dust 
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Fig. 27. Histograms of the 353 GHz dust unit conversion / colour 
correction maps for the Full survey (a.), Nominal survey (b.), 
survey 1 (c.) and survey 2 (d.). The plots include histograms 
based on the full sky as well as the LMC, CHA, and NEP re- 
gions. 



emission is much stronger with higher frequencies, and is thus 
less dominant at 143 and 100 GHz. 

A dust colour correction coefficient originating from the CO 
extraction (Planck Collaboration XIII 2013) is also shown for 
the 217 and 353 GHz bands. 

The 217 and 353 GHz portions of Fig. 28 also contain data 
points from the CO extraction (Planck Collaboration XIII 2013). 
The horizontal bars on some of these points indicate that the 
coefficient was derived for the PSB A/B pair of detectors, not the 
detectors individually. For the 217 and 353 GHz coefficients, the 
bandpass and CO values are in excellent agreement, even where 
the dust coefficients from this study appear to diverge from the 
bandpass values slightly. The exact causes of these variations 
remains under study for a future data release. 

The uncertainty in the figures for the sky-based dust coeffi- 
cients is based on an absolute calibration uncertainty of 0.5 % 
for each band. The bandpass based coefficient uncertainties are 
based upon the quoted uncertainties of and p of 1.9 K and 
0.25, respectively, as well as the spectral uncertainty associated 
with each spectral response profile. The CO coefficient uncer- 
tainties are based upon a 1 % relative calibration uncertainty. 

4.2.3. S-Z Bandpass Verification 

The Sunyaev-Zeldovich (SZ) eff'ect (Sunyaev & Zeldovich 
1980) has a characteristic signature in the millimetre/ submil- 
limetre domain. With a changing sign at 217 GHz, it is quite 
diff'erent from the CMB anisotropy spectrum and other power- 
law spectra typical of these frequencies. Consequently, any 
bandpass leakage should be clearly identifiable in SZ spec- 
tra. To isolate the SZ signature from other foregrounds and 
the CMB anisotropics themselves, the 20 brightest clusters in 
the Planck catalog (Planck Collaboration XXVIII 2013, Planck 
Collaboration XXIX 2013, and Planck Collaboration XXI 2013) 
have been selected to provide lines-of- sight over which to in- 
vestigate spectral mismatch leakage. The lines-of-sight selected 
are A2319, A3266, RXC J1638.2-6420, A2219, A2142, Coma, 
A366, A2255, A2029, A3186, A2218, A3158, A85, A3827, 
A697, A1795, A644, A2204, A3628, and A3888. For each of 
the above sources, the integrated flux in individual HFI chan- 
nel maps is measured (calibrated in thermodynamic temper- 
ature KcMB units) at the position of each cluster. Using the 
nominal bandpass conversion coefficients for each of the ex- 
pected components, the integrated 7sz for each cluster is de- 
duced from the residuals using a statistical analysis. Thus, 
an estimate of the SZ unit conversion from the Fsz Compton 
parameter to the <ircMB difl'erential CMB temperature may be 
obtained for each cluster, and each frequency; without invoking 
the bandpass-based SZ unit conversion coefficients themselves 
(Eq. 33, Table 6). This analysis provides a consistency check of 
the bandpass-based SZ coefficients. 

Two methods have been used to perform this analysis. The 
first allows only a dust template removal, fixed on the 857 GHz 
channel, with coefficients deduced outside each cluster. This first 
method does not make assumptions about the colour of the back- 
ground. The background removal was accomplished via subtrac- 
tion of an average value taken on an annulus surrounding the 
cluster, i.e., aperture photometry. This method has a low signal- 
to-noise ratio because it is dominated by the CMB anisotropy 
residuals. The second method uses the MILCA algorithm (Hurler 
et al. 2010) to subtract a component with a CMB spectrum at all 
frequencies. The average of these coefficients per frequency is 
shown in Table 8 for the two methods and compared with the ex- 
pected coefficients based on the nominal band-average transmis- 
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Table 8. Various jsz to dTQy[^ unit conversion factors. 



Band 


^cBP 


Ap. Phot 


^cMILCA 


[GHz] 


[jsz/Kcmb] 


[jsz/Kcmb] 


[jsz/Kcmr] 


100 


-4.030 + 0.018 


-4.7 + 0.6 


-4.6 +0.2 


143 


-2.78 +0.04 


-3.9 + 0.6 


-3.00 + 0.10 


217 


0.19 +0.05 


-1.3 + 0.6 


-0.10 + 0.10 


353 


6.21 +0.11 


4.3 + 0.8 


5.6 +0.2 


545 


14.46 +0.07 


17.4 + 6 


17 +2 



sion spectra. Although the signal-to-noise ratio is not larger than 
typically 10 for the first method and 30 for the second method, 
the agreement is remarkable. Thus, the SZ cluster data provides 
a pseudo-quantitative verification of the accuracy of the HFI 
detector spectra as gross spectral leakages would demonstrate 
themselves through this analysis. 



^ 5. Conclusions 



Fig. 28. Comparison of the sky and bandpass based dust unit 
conversion / colour correction coefficients for the HFI 100 GHz 
-353 GHz spectral bands (a. - d.). 



The spectral response for the HFI detectors has been presented. 
The derivation of the HFI frequency channel band-average spec- 
tra has been presented, including photometric/noise/sky cover- 
age constituent scaling coefficients. The scaling coefficients for 
the band-average spectra were compared with those resulting 
from the various individual HFI surveys, where it was demon- 
strated that the individual survey values converged to the full 
survey average values used in the derivation of the spectral re- 
sponse data products. Unit conversion and colour correction co- 
efficient relations have been derived. The corresponding coef- 
ficients, and uncertainties, are to be made available within the 
Planck data archive (Planck Collaboration ES 2013). The accu- 
racy of the HFI spectral response has been verified using both 
ground-based component level data and, importantly, HFI flight 
data. The defining requirement indicated in the original HFI cal- 
ibration plan is knowledge of the spectral transmission of the 
individual detectors with uncertainties below 3 % for the low 
frequency channels (100, 143, and 217 GHz) and below 1 % for 
the high frequency channels (353, 545, and 857 GHz); all with 
a spectral resolution of better than 3 GHz. The spectral resolu- 
tion requirement has been exceeded by more than a factor of 
five. It is possible to degrade the spectral resolution requirement 
to gain an improvement in the S/N, which allows the desired 
spectral transmission accuracy to be achieved in most cases (the 
absolute uncertainty for the 100 GHz detectors is high due to the 
high noise levels of the reference bolometer at these frequen- 
cies, but the relative uncertainty even for the 100 GHz detectors 
remains low). Estimates of the out-of-band transmission profiles 
have been incorporated into the bandpass data products. Out- 
of-band signal attenuation is demonstrated to be better than 10^ 
through HFI observations of the diff'use zodiacal cloud Planck 
Collaboration XIV (2013). Good agreement has been demon- 
strated between sky and bandpass dust colour correction coeflft- 
cients, as well as SZ unit conversion coeflftcients. The sky and 
bandpass based CO coefficients have shown a correlation, yet a 
full comparison is unconstrained at this time due to diff'erences 
between the two approaches. The CO comparison is a work in 
progress, where improvements are expected as Planck polariza- 
tion data analysis progresses. Software tools to calculate Planck 
unit conversion and colour correction coefficients have been pro- 
vided, along with tables of coefficient values for individual de- 
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tectors, and band-average spectra, in Planck Collaboration ES 
(2013). 
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